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PROGRESS . REPORT 

This report summarizes research carried out at Harvard 
University under NASA grant NGL 22-007-247 during the report 
periods 1 August 1975 through 31 January 1976 and 1 February 1976 
through 31 July 1976. During this period eight papers have been 
published or accepted for publication; five of these were written 
during the previous reporting period and three papers report new 
work completed this year. In addition, nine talks on lunar re- 
search were given at professional meetings; all but one were 
piiblished in abstract form. 

During this report year our research has followed generally 
the lines set forth in last year's report and proposal. The cooling 
rate studies described in the earlier report have now been applied 
to the Apollo 12 mare basalt suite in an attempt to deduce the geometry 
of the* flow or flows from which they came. This involved using the 
experimental data on crystallization history and liquid compositions 
to calculate liquid viscosities and densities as a function of time, 
then making model calculations of crystal settling during cooling 
of a magma body, and comparing the results obtained with observations 
on the Apollo 12 samples. We conclude that the samples studied ' 
came from near the bottom of a rather thick flow. This work is 
described in paper 7 (preprint attached) . 

Study of titaniferous basalt samples 74275 and 70215 revealed 
a significant degree of heterogeneity within each of these rocks 
in terms of- major - as well as minor - element composition. The 
reason for this heterogeneity is not yet fully understood, but it 
must be taken into account in discussions of the origin of. these 
rocks. We have shown that the observed heterogeneity is sufficient 
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to accoxmt for the apparent conflicts between the experimental 
crystallization data on these two samples^ and that the existence 
of this heterogeneity severely limits our ability to determine 
tHe depth-of-origin from phase-eguilibriuin- data. See paper 6 
(reprint attached) for details. 

Our efforts to codify and summarize the very large amount of 
chemical data on experimental run products continue. The relation- 
ship between Fe and Mg in olivine and cbexisting liquid was 
described last year, generally confirming the earlier results of 
Roedder and Emslie who used terrestrial materials.. This year 
John Longhi completed a lengthy study of Fe and Mg in plagioclase 
and the related problem of stoichiometry in Ixmar plagioclase. 

The essence of this work with lunar applications is reported in 
paper 8 and abstract 6 (preprints enclosed) . Further details are 
in Longhi *s Ph.D. thesis. 

i 

During the course of the year. Walker designed and built a 
molybdenum strip furnace for the melting of rock samples to 
facilitate XRF and microprobe analysis for major and- minor element 
bulk compositions. Although this apparatus was not intended to 
be used with l^xnar samples, some interesting lunar applications 
have been found. For examplej it was easy to show that calcic 
pi-agioclase readily floats in a synthetic melt with a composition 
matching that of the liquid from which lunar crustal plagioclase 
is believed to have crystallized. This puts to rest a question 
that has been hotly debated (see abstract 9, attached) , 

Furthermore, with the strip furnace it may be possible to 
obtain heats of fusion for lunar rocks or their analogs. Such 
data is needed for many purposes, and is almost totally lacking 
for terrestrial as well as lunar rocks. Work is in progress to 
develop and prove out the technique. 
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JohJi Longlii received liis Ph.D. degree iii June 1976 and is now 
a Postdoctoral Research Fellow at MIT'. Edward Stolper has received 
his M.Phil. degree from the University of Edoiiburgh for his thesis 
on the experimental petrology and origin of basaltic achondrites 
and will return to Harvard in September to rejoin the lunar group 
here. 
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PUBLICATIONS 

{1 August 1975 - 31 July 1976) 

Articles in Boolcs or Scientific Journals 

1. Lunar igneous rocks and the nature of the Itinar interior. 

J.F. Hays and D. Walker. Proc. Sbviet-American Conf. on 
Cositiochemistry of the Moon and PlanetSj Moscow. The Lunar 
Science Institute (in press). 

2. Origin of titaniferous lunar hasaits, D. Walker. J. Longhi, 

T.L. Grove, E.M. Stolper, and J.F. Hays. Geochim, Cosmochim. 

Acta, 3^ 1219-1235, 1975. 

3. Differentiation of a very thick magma body with implications 
for the source regions of mare hasaits. D, Walker, J. Longhi, 
and J.F. Hays. Proc- Sixth Lunar Sci. Conf., Geochim. 

Cosmochim. Acta, Suppl. 6, 1103-1120, 1975. 

4. Crystallization history of lunar pier itlc basalt sample 12002: 
phase-equilibria and cooling, rate studies. D. Walker, R.J. 
Kirkpatrick, J. Longhi, and J.F. Hays. Geol. Soc. Amer. Bull., 

87. 646-656, 1976. 

5. Direct determination of the quart z-coesite transition by in-situ 

x-ray measurements: a discussion. J.F. Hays, Contrib. Mineral- 

Petrol., ^ 61-63, 1975. 

6. Heterogeneity in titaniferous lunar basalts. D. Walker, J. 

Longhi, and J.F. Hays. Earth Planet. Sci. Letters, 30, 27-36, 1976. 
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7. Differentiation of an Apollo 12 picfite magma. D, Walkerj^ 

J. Longhi, R.J. Kirkpatrick, and J.F. Hays. Proc. Seventh 
Lunar Sci. Conf., Geochim. Cosmochira. Acta, Suppl. 7 (in press). 

I 

B. Fe and Mg in plagioclase. J. Longhi, D. Walker, and J.F. Hays. 
Proc. Seventh Lunar Sci- Conf., Geochim. Cosmochim. Acta. 

Suppl. 7 (in press). 
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Abstracts a nd Tallcs given at Meetings 

1. Acceptance speech for F.W. Clarke Medal. D. Walker. Geochlm. 
Cosmochim, Acta, 40_, 567-569, 1976. 

2. Heterogeneity in titaniferous lunar basalts. D. Walker, 

J. Longhi, and J.P. Hays. Origins of Mare Basalts. The Lunar 
Science Institute, Houston, Texas, 169-173, 1975. 

3. Experimental petrology of mare basalts: an introduction. 

J.F. Hays. Given orally at .Mare Basalt Conference, The 
Lunar Science Institute, Houston, Texas, November 1975. 

4. Differentiation of a thin magma body. D. Walker, J. Longhi, 
and J.F. Hays. Lunar Sci ence VII. 904-906, The Lunar Science 
Institute, Houston, Texas, 1976. 

5. Fe, Mg, and silica in lunar plagioclase. J. Longhi, D. Walker, and 
J.F. Hays. Lunar Sci ence VII. 501-503, The Lunar Science 
Institute, Houston, Texas, 1976. 


6. Solid solution and phase equilibria on the join anorthite- 

silica. J. Longhi and J.F. Hays. Trans. Am. Geophys. Un. (E®S), 
57., 340, 1976. 


7. Olivine nucleation in lunar basaltic melts'’’ in iron capsules. 

D. Walker, R.j. Kirkpatrick, J. Longhi, and J.F. Hays. Trans. 
Am. Geophys. Un. (E#S), ^ 356, 1976. 

8. Crystallization of basalts in the laboratory. J.F. Hays. Geol. 
Soc. Amer. 1976 Ann. Meeting' Rocky Mountain Section, Abstr. 
with Prog. 8 ^ 590, 1976. 
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9. Feldspar flotation and Itmar crust foirmation* D. Walker and 
J.F. Hays. (submitted for presentation at Geol. Soc, Amer. 
1976 Ann. Meeting) . 
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HETEROGENEITY IN TITANIFEROUS LUNAR BASALTS 
DAVID WALKER, JOHN LONGHI and JAMES F. HAYS 

Hoffman Laboratory; Center for Earth and Planetary Physics. Harvard University, 
Cambridge, Mass. (USA) 

Received December 12, 1975 
Revised version received January 28, 1976 


Small, but real, chemKal differences exist between subsamples of finc-greined, qucnch-textured titaniferous lunar 
basalts. The existence of different textural domains with diffexent chemistries is thoiigbt to account for most of this 
variation. In addition to the textural domains, 74275 has a population of a few percent of F 071.80 olivine "mega- , 
ctysts”, as well as dnnite fragments of Fo 80 l 84 - These materials are thought to be extraneous and to compromise the 
primary nature of 74275. 

Recognition of the small chemical variations present may aid in understanding some discrepancies in the expai- 
mcntal oetiology iiteiature. However, these small variations have a distressing petrogcnetic significance since they 
severely limit resolution in recognizing the number and depth of origin of primary magmas. 


1. Introduction 

Titaniferous lunar basalts returned from the Apollo 
1 1 and 1 7 landing sites have proved to be a continuing 
subject of controversy in lunar petrogenesis. Originally 
the Apollo 1 1 basalts were interpreted on the basis of 
experimental petrology either as extreme residues of 
fractional crystailiaation near the lunar surface [1] or 
as incipient melts of the primitive lunar interior [2]. An 
alternative interpretation based on trace element 
studies [3] was that the titaniferous basalts represented 
the melting products of evolved regions of the lunar 
interior. The study of fine-grained, igneous textured 
titaniferous basalts from Apollo 17, which could not 
be explained as the residue of near-surface fractional 
crystailiaation, produced experimental evidence that 
the parental titaniferous basalts could be produced at 
relatively shallow depths (100-i 50 km) in the moon 
by melting of an flmenite-bearing source region [4,5].' 
This evidence supported the interpretation of melting 
of cumulates and was consistent with the idea that the 
moon had been melted to a depth of at least a few 
hundred kilometers very early in its history and had 
differentiated into a series of floating and sinking 
cumulate layers [6—9]. In this model the mare basalts 


are derived as the melting product of the sunken mafic 
cumulates. The principle difference between low- and 
high-Ti basalts is that the high-Ti basalts are derived 
from shallower depths by melting of the last formed, 
most differentiated cumulates whereas the low-Ti ba- 
salts have their source regions in the deeper, earlier 
cumulates. 

This interpretation of the high-Ti mare basalts is 
based on finding ilmenite, olivine, and clinopyroxene 
on the liquidus of 70215, a fine-grained igneous sam- 
ple, at between 5 and 7.5 kbar. This interpretation has 
been challenged [10] on the basis that 70215 is not the 
best representative of parental high-Ti basalt. 74275 
with its higher Mg value (100 X Mg/Mg + Fe molar) 
and more olivine-normative character was alternatively 
proposed [iO] as a parental high-Ti magma. Smee 
74275 had. a higher pressure of olivine and pyroxene 
cosaturation (12 kbar) and since ilmenite is not also a 
saturating phase, the incentive for generating high- and 
low-Ti basalts from different source regions at different 
depths was greatly reduced. 

It is the purpose of this paper to show that 70215 
and 74275 are not so dissimilar as the previous studies 
suggest [4,10]; tlat the small differences result from 
sample heterogeneity coupled with minor olivine accu- 
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rmilation;an<i that, while recognition of the differences 
can account for supposed experimental discrepancies, 
the differences are not sufficient to require drastic 
changes to existing models for the petrogenesis of high- 
Ti basdts 


2. Textural domains and chemical variations 

Rocks 70215 and 74275 are finergramed, rocks 
with petrographic characteristics indicating rapid c^s- 
tallizalion from a silicate liquid. These characteristics 
include skeletal microphenocrysts, feathery to spheru- 
litic crystal intergrowths, elaborate chemical zonation 
in minerals, and preservation of metastable armalcolite 


and mesostasis glass in the solidification product. It 
was therefore assumed that the chemistries of both 
70215 and 74275 represented the chemistries of liq- 
uids which existed at the surface of the moon and 
petrogenetic discussion proceeded on that assumption. 

Subsequent analytical and experimental work have 
cast doubt on this assumption, inasmuch as significant 
chemical heterogeneities appear to exist among differ- 
ent subsamples of each rock [28,24], Figs. 1 and 2 
display the published analyses of subsamples of 70215 
and 74275 [10-18]. These projections have proved to 
be useful in discussing the crystallization behavior of 
titaniferous lunar basalts [5]. It is dear that a range of 
compositions is present in both 70215 and 74275 and 
that the ranges are about the same size in each case. 



Fjg. 1. jProjection astern for titaniferous basalt compositions used in fig. 2 of Walker et al. [5]. Underlined minerals are the liq- 
uidus phases in that part of composition space. Shaded blob covers compositions of low-K, high-Ti basalts of Apollo 1 1 and 1 7. 
Enlarged triangle shows reported compoations of 70215 and 74275 [10-18]. Fuzzy and giainy are compositions of textural do- 
mains from Table 1. C and H are compositions of 70215 studied at Canberra [10} and Harvard [4], Rosette of arrows indicates the 
magnitude and direction of displacement of compositions corresponding to the analytical precision of each oxide. 
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Fig. 2. Projection system for tttaniferous basalt compositions from FmTiOs used in fig. 3 of Walker et al, {5j. Underlined minerals 
are tlie Uquidus pliases in that part of composition space. Shaded blob covers compositions of low-K, high-Ti basalts of Apollo 1 1 
and 17. Enlarged triangle shows reported compositions of 70215 and 74275 {iO-lSj. Fuzzy and grainy are compositions of 
textural domains from Tabic 1. Rosette of arrows indicates the magnitude and direction of displacement of compositions corre- 
sponding to the analytical precision of each oxide. 


The variations are actually rather srnall, however, they 
are significantly above the analytical precision (about ' 
a factor of 3). This is especially clear since two sub- 
samples of 70215 were analyzed in the sante laboratory 
by.the same techiiique'[16,17j. It is concluded that 
real chemicai differences exist between subsamples of 
70215 and 74275. Furthermore the differences cannot 
be the result of simple crystal settling or local fraction- 
al crystallization. The chemical variation does not con- 
'form to simple olivine or armaicolite fractionation 
curves. In Fig. 1 the cosaturation liquidus curve of 
olivine with armaicolite [4] cuts across the field of 
70215, If the range in 70215 were caused by local 
fractional crystallization, the variation should follow 


tliis liquidus boundary and not cross it. 

How, then, do these chemical differences arise? 

Fig. 3A shows section 70215,147 to be texturally 
heterogeneous. Fig. 3B shows in an enlargement of 
70215,147 that the textural domains differ principally 
in crystal morphology. The more translucent domain 
is characterized by better crystallized, lightly pink, 
granular pyroxene whereas tlie darker domain is char- 
acterized by more poorly crystallized spherulitic inter- 
growths of fuzzy-appearing pyroxene with the other 
minerals. These domains are dubbed “fuzzy” and 
“grainy” for convenient reference. Similar domains are 
found in 74275 as well as in other high-Ti Apollo 17 
basalts (e.g. 71569, 72155, 72135). Fig. 3C shows 
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Fig. 3. A. 70215.147; transmitted light, note heterogeneous textural domains. 13. 70215,147; enlargement of A. C. 70215,156; 
transmitted light. D. 70215,156; enlargement of C, note skeletal olivine microphenocrysts. 


TABLE 1 


Textural domain and olivine compositions 



70215,147 

70215,147 

74275,97 

Olivines ♦* 



grainy ♦ 

fuzzy • 

(a) 

(b) 

(c) 

SiOz 

38.0 

37.2 

38.8 

38.8 

37.7 

TiOz 

12.8 

13.3 

0.13 

0.13 

0.18 

AlzOj 

9.2 

8.6 




CrzOa 



0.31 

0.37 

0.28 

FeO 

18.7 

19.9 

15.6 

17.6 

22.3 

MgO 

8.0 

8.7 

44.0 

41.8 

38.4 

MnO 



0.20 

0.24 

0.28 

CaO 

10.8 

10.6 

0.26 

0.25 

0.32 

NiO 



0.06 

0.02 

0.05 

ZnO 



0.02 

0.17 

0.21 


97.5 

98.3 

99.58 

99.38 

99.72 


Electron microprobe defocused beam partial analyses. Oxide percentages in each domain, after normal matrix correction proce- 
dures, were normalized to the average of published 70215 oxide values to correct for abnormal matrix effects inherent in defocused 


beam techniques applied to analysis of thin sections. Values used for 70215 average, grainy, and fuzzy are; SiOz, 37.61, 43.2, 42.3; 
TiOz, 13.05, 11.1, 11.5; AlzOz, 8.88, 12.7, 11.9; FeO, 19.34, 17.6, 18.7; MgO, 8.35, 7.6, 8.2;CaO, 10.69, 11.1, 10.9. 

(a) Core of dunite xenolith in Fig. 3H. (b) Core of large euhedral phenocryst. (c) Microphenocryst in groundmass. 
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Fig. 3 (continued). E. 74275,94; transmitted light, crossed polarizers. Olivine “megacryst”. F. 74275,94; transmitted light, 
crossed polarizers. Olivine xenolith. G. 74275,97; transmitted light. Large euhedral olivine crystals. H. 74275,97; transmitted 
light, crossed polarizers. Dunite xenolith with annealed granular texture. 


that these domains are not small in comparison with a 
thin section (and therefore a subsample for analytical 
or experimental work) since 70215,156 is predomi- 
nantly “grainy” domain texture. Defocused beam 
microprobe analysis was done on different domains in 
70215,147 (see Table 1). The results are indicated in 
Figs. 1 and 2. The domains have measurably different 
chemistry and vary in the same sense as the range of 
compositions reported for the subsamples. Thus the 
inclusion of variable proportions of these domains in 
the subsamples appears to be capable of causing the 
compositional variation among the subsamples of a 
given rock. Although the reported subsamples show a 
slightly larger range in chemistry, it is thought that the 


different textural domains occurring in different pro- 
portions explains much of the observed chemical varia- 
tion between subsamples. 


3. Origin of textural domains and implications 

The origin of the textural domains is unknown. It 
is possible that the liquid was heterogeneous and that 
the different textures result from crystallization of 
parts of the liquid with different chemistry. This is 
considered to be extremely unlikely because textural 
differences appear to be much more strongly control- 
led by cooling history than by chemistry (compare 
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[19,20]) and the chemical differences are small. Fur- 
thermore, considering the ease with which these low- 
viscosity melts homogenize during laboratory fusion, 
it seems unreasonable to expect that heterogeneities 
could be maintained in the liquid for any appreciable 
time prior to solidification. This suggests that the 
heterogeneities are introduced during the solidifica- 
tion of a homogeneous liquid and requires that the 
material transfer rate within the solidifymg magma 
exceeded the solidification rate. The need for some 
process to accomplish this objective has been antici- 
pated in previous discussions of lunar basalt petrogene- 
sis [21]. Briefly, the difference between the composi- 
tion of plagiodase-saturated residual liquids and hand 
specimen compositions has been argued to result 
either from microphenocryst accumulation in the hand 
specimen (for wliich evidence is poor) or from some 
other unrecognized mass transfer process. In the pres- 
ent case it has been suggested that the granular areas 
may represent globules of solidified material-which 
sank into material of slightly difTerent composition 
which ultimately crystallized as the fuzzy textural do- 
main [22]. The mass transfer distance might be char- 
acterized as “large” compared to a hand specimen. 

Alternatively, it may be that the development of 
grainy and fuzzy areas is due to nudeation and growth 
effects under conditions of rapid cooling. It has been 
proposed [4] that early formed crystals nucleated and 
began to grow rapidly in isolated clusters with release 
of latent heat which locally reduced the codling rate. 
This allowed pyroxene to nucleate before plagiodase 
and to coarsen, and permitted some olivine resorption. 
Between these dusters, nudeation of pyroxene was 
delayed until greater undercoolings were reached and 
the nudeation rate was much higher. As a result, when 
nudeation did take place, pyroxene and plagiodase 
together with ilmenite formed very fine radiate inter- 
growths, which now have a fuzzy appearance, and 
there was little opportunity for olivine resorption. 
Compositional variation among the domains results 
from the differential incorporation and exclusion of 
material as the clusters grew. Presumably there were 
compositional gradients in the melt between the 
clusters. In this model the mass transfer distance 
might be characterized as “small” compared to a hand 
specimen since material exchange occurs only on the 
scale of the domains. 

Whatever the origin of these textural domains, 


their demonstrated existence and chemical variation 
has several important implications. First it is fruitless 
to argue in detail about differences in experimental 
crystallization sequences and liquidus temperatures of 
allocated subiamples since these subsamples quite ob- 
viously are heterogeneous. In particular in Fig. 1, note 
that 70215 compositions cross the olivine-armaicolite 
liquidus boundary. Some subsamples should crystallize 
olivine first and some should ctystailize armalcolife 
first, in iron capsules in a sealed tube. Composition H 
studied at Harvard [4] crystallizes armalcofite first and 
C studied at Canberra [10] crystallizes olivine, armal- 
cohte (and spinel), at the liquidus. A synthetic compo- 
sition was studied at Stony Brook [23] which crj'Stai- 
lized both minerals simultaneously, althougli the re- 
potted mix did not lie exactly on the olivine-armalco- 
lite curve shown here. It is thought that allocation 
heterogeneity may explain some of these experimen- 
tal discrepancies. It should also be noted here that the 
large differences implied to exist between 74275 and 
70215 [iO] may not be so significant as thought Such 
conclusions are strongly dependent on particular sub- 
samples studies. The Mg number appears to be the 
major difference between the two samples, as will be 
discussed below. 

Clearly subsampies of 70215 or 74275 may not re- 
present liquid compositions as such, although they 
may have averages which were once liquids. Pieces of 
the same rock can be unrelatable to one another by 
normal crystal-liquid processes and might inspire a 
conclusion of a separate petrogenesis. Such a conclu- • 
sioa would be spurious. A similar state of affairs may 
also apply to some subsamples on which trace elements 
have been measured [24], On the scale of investigation 
implied in sample allocation there is a danger of infer- 
ring more primary magmas than actually existed. 

A final implication may be noted about the magni- 
tude of chemical variation observed in the textural 
domains. They were small and we conclude that they 
result from mass transfer induced during solidification 
of a liquid. The scale of the mass transfer is not resolved 
and it could conceivably be large. Evan if this is so, it 
must be noted that the variations documented are 
quite insufficient to explain the difference between 
hand specimens and cotectic liquid chemistry. Some 
additional mass transfer process would be required to 
explain the difference. We prefer to regard the cotectic 
liquids sampled (e.g. 75035) as products of fractional 
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crystallization of parental hand specimen composi- 
tions such as 70215 or 74275 since there is a continuous 
suite of hand specimen samples whose compositions 
follow a liquid line of descent from such parental com- 
positions [4,5]. The compositional variation of the 
whole suite would require some rather improbable 
coincidences if they were to be interpreted as the prod- 
ucts of crystal accumulation. 

Note in Fig. 2 that the composition of the low-K, 
high-Ti basalt suite shows a non-ti'near pattern. It has 
been documented [4,5,10] that this “bent’’ pattern 
conforms to a liquid line of descent. It is difficult to 
see'how this “bent” pattern could be produced during 
crystal accumulation in plagioclase-saturated liquids. 
Furthermore the accumulation model involving plagio- 
clase-saturated liquids fails because the armalcolites 
actually found in the samples are not stable in pyrox- 
ene-plagioclase-saturated liquids [4,5,10,23]. We do not 
wish to preclude the possibility that some other mass 
transfer process may be operative, but crystal accumu- 
lation is highly improbable and whatever process leads 
to textural domain formation is quantitatively insuffi- 
cient. 


4. 70215 versus 74275 . 

An additional feature of the petrography of 74275 
needs discussion. Both 74275 and 70215 have skeletal 
microphenocry'sts which occasionally are as large as 
300 nm. See Fig. 3D. They arc generally Fo,, or less 
magnesian. In contrast to any other high-Ti basalt thin 
sections we have seen, thin sections of 74275 have a 
distinctive population of olivine “megacrysts” wliich 
reach 3 mm in size (Fig. 3E). Their very irregular 
morphology suggests that they are not phenocrysts 
trapped in a chilled liquid. Furthermore this popula- 
tion has occasional “megacrysts” in which more than 
one crystal is present very much as if it were a sm^ 
xenolithic fragment of an olivine-rich rock with an- 
nealed granular texture (Fig. 3F). It is thought, on 
textural grounds, that these “megacrysts” do riot re- 
present phenocrysts grown from a 74275 melt but 
rather are crystal fragments included in 74275. Hiis - 
population comprises only a few percent of 74275 
and compositionally is relatively magnesian Fovf-so- 
Section 74275,97 has some special features in addi- 
tion to these anhedral “megacrysts” [29]. It contains 


a few large (up to 1.5 mm) euhedral olivine crystals 
with cores as magnesian as Fogi (Fig. 3G). Also pres- 
ent is a 1 -mm polycrystaOine olivine clast with an- 
nealed granular texture (Fig. 3H). Crystals in the cen- 
ter of the clast are Fog 4 whereas those on the edge are 
as iron-rich as F 074 . This clast is interpreted as a xeno- 
lith of dunite [29] which has suffered some Fe-Mg 
exchange with the melt along its margins. 

74275 has been claimed to be significantly more 
primitive (higher Mg-value and more normative olivine) 
than 70215 [10]. We have noted above that the magni- 
tude of the presumed difference between 70215 and 
74275 is quite dependent on which subsamples are 
examined and that in Figs. 1 and 2 there is quite a 
close approach of some of the subsamples of the two. 
Unfortunately these projections suppress Fe— Mg varia- 
tion so this must be displayed separately in Fig. 4. The 
range of published compositions does not overlap here, 
thus.showing that 74275 and 70215 are indeed dis- 
tinguishable on this basis. Also shown m Fig. 4 are 
curves for oiivme-liquid distribution coefficients of 
Fe-Mg forifo^ 0-26 and 0.28, which are thought to 
be relevant to the crystallization of these titaniferous 



Fig. 4. Molar Fe/Fe + Mg of reported analyses of 70215 and 
74275,110-18}. Also plotted are Fc/Fe + Mg values of most 
magnesian olivines found in each sample, h’or 74275 this is 
taken to include the iaige euhedrai olivines with phenoeryst 
morphology (Fig. 3G) but does not include the F 084 found in 
the dunite xenolith (Fig. 3H). Also shown are two oiivlne/liq- 
uid distribution coefficient curves for Fe/Mg f25,10}. 
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lunar liquids [25,10]. It can be seen that within the 
precision of the data, the most magnesian olivines pres- 
ent in 70215 could have grown directly from liquids 
of 70215 composition. This does not appear to be 
true for 74275. Clearly the Fos 4 present in the dunite 
xenohth could not have grown from liquids of 74275 
composition. The Fosi shown in Fig. 4 is for cores of 
large (millimeter size) euhedral olivines which have 
phehocryst morphology, hdwever, it can be seen that 
even these are also too magnesian to be interpreted as 
the direct crystallization product of liquids of 74275 
composition. Furthermore we questioned above on 
textural grounds whether the magnesian anhedral mega- 
crysts in 74275' really do represent crystallization prod- 
ucts of a 74275 liquid. We point out that exclusion of 
only 4% of such “megacrysts” from 74275 produces a 
composition with Mg value in the 70215 range. That 
is to say if 4% olivine megacrysts are added to a liquid 
with Mg value of 70215, the resultant mixture would 
have the higher Mg value and more normative olivine 
characteristic of 74275. This would require that the 
megacrysts and other peculiar olivines be xenocrysts 
from an unknown source, since olivines of this compo- 
sition should not have grown from a liquid with 70215 
Mg value. The xenocryst hypothesis appears to con- 
forhi to some of the petrographic features of these 
peculiar olivines noted above. 

We are not trying to prove that [olivine + 70215 
= 74275] actually occurred since trace element char- 
acteristics seem to rule out such a direct relationship 
[23]. However, in view of the major element hetero- 
geneities between subsamples, this hypothesis may not 
be impossible. Evaluation may critically depend on the 
choice of subsamples. 


5. Petrogenetic discussion • 

A substantial part of the incentive for interpreting 
high-Ti basalt parental liquids (such as 70215) as the 
fusion product of Ti-enriched cumulates at relatively 
shallow depths within the moon was the discovery 
that 70215 was multiply saturated with olivine, pyrox- 
ene, and iimenite at the liquid us at between 5 and 7.5 
kbar [4]. This observation has been challenged [10] on 
the basis of experimental work on a different* subsam- 
ple of 70215. This work did not show multiple satura- 
tion with all three phases at '5 or 8 kbar although ilme- 
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Fig. 5. Experimental results. on 70215,69 as a function of 
pressure and temperature. Ail experiments in high-purity ircm 
capsules. 70215,69 is saturated with olivine, pyroxene, and 
ilmemte between 5 and 7.5 kbar. Low-pxessuie crystallization 
from [4]. ■ 


nite was dose to olivine and pyroxene at S kbar. Even ■ 
if the subsaraples had the same chemistry, which they 
do not, it seems imprudent to us to challenge the ori- 
ginal observation without making runs in the suggested 
pressure range between 5 and 7.5 kbar. We have per- 
formed new experiments in high-purity iron capsules 
to check our previous results in molybdenum and 
graphite [4] (see Fig. 5). At 1205°C and 7.5 kbar, 
70215;69 is saturated with pyroxene and ilmemte less 
than 10®C from the iiquidus; at 1 190“C and 5 kbar, 
again less than 1 0°C front the Iiquidus, it is saturated 
with olivine alone. Spinel does not occur in this multi- 
ple saturation assemblage when Fe capsules are used. 
We feel that this confirms our earlier work showing 
that 70215,69 is saturated with olivine, pyroxene, and 
iimenite at ,6-7 kbar. Therefore, if 70215,69 is repre- 
sentative of a melt composition generated in the lunar 
interior, it may be interpreted as the melting product 
of an iimenite-bearing source region at 100—150 km. 

Since heterogeneities do exist between subsamples 
it is necessary to inquire how sensitive the multiple- 
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FeO , MgO 
71.8 40.3 


Fig. 6. (FeO/71.8 + MgO/40.3) is plotted versus tlie Hquidus 
cosaturation pressure of olivine and pyroxene for lunar basalt 
compositions. Data taken from references 4, 5, 1 0, 20, 26 and 
27. Different lines wuh the same slope apply to high and iow- 
Ti compositions. “Error bars” on 70215 and 74275 indicate 
the range in repoiied analyses. Tliis variation causes an uncer- 
tainly of ~75 km in estimating the depth to olivine + pyrox- 
ene cosaturation for 70215. 


saturation pressure (and hence depth-of-origin estimate) 
is to these small chemical variations. The pressure of 
cosaturation with pyroxene and olivine shows a regular 
dependence on the parameter (FeO/71.8 + MgO/40.3) 
which is a measure of olivine abundance. Fig. 6 shows 
that different lines having the same slope apply to 
high-Ti and low-Ti lunar compositions. The variations 
in composition discussed above for 70215 and 74275 
are indicated by the “error bars” in Fig. 6. Clearly the 
depth estimate is quite sensitive to these small chemi- 
cal variations. The estimate for 70215 may jvaiy by as 
much as 75 km based on the differences between sub- 
samples alone! It is to be noted that if 4% “megacrysts” 
are removed from 74275 then the inferred depth of 
origin becomes the same as that of 70215, i.e. 100—1 50 
km. 


6. Summary 

Perhaps it is pointless to argie in detail about which 
subsamples represent primary magmas, at what depth 


they are generated, and how many are necessary. Cer- 
tainly the factor of sample heterogeneity must be 
neutralized in some way before further progress can be 
made. How this is to be done with the present data 
baseds unclear. However, at the present state of our 
investigations we can probably rule out “parental” 
samples such as 70215 as being products of crystal ac- 
cumulation into.cotectic liquids. If the “parental” 
compositions are in fact “primary” then the shallow- 
melting-of-Ti-rich-cumulates scenario follows. Rejec- 
tion of this scenario, on the basis of results from 
74275, is premature. 
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Abstract 

Dlabributlon coofficlonta for Fe and between plagioclaee 
and basaltic liquids have been measured for lunar, terrestrial 
and synthetic systems. Observed substitution in plagioclase 
is consistent with a Ca(Fa,Mg)Si20g component, and Mg is 
incorporated into plagioclase more readily than ferrous iron. 

for lunar basalts and synthetic analogs. As 
oxygen fugacity increases, ferric iron becomes significant and 
the bulk exchange coefficient Kp'j^ increases to an average . 
of 1,4 for. terrestrial basalts. The variation of with 

oxygen fugacity provides a crude oxygen barometer. Lower con- 
centrations of Fe and Mg in plagioclase from plutonic rocks 
suggest that orystalllaation kinetics significantly effect tho 
observed distributions. 

The new data imply tbat<- 1) lunar anorthosites. crystallized 
from relatively iron-rich liquids, 2) Fe/Mg zoning in plagioclase 
is an Indicator of crystallization processes in basaltic 
3) Fe/Mg heterogeneity in plagioclase 1s a clue to the recognition 
of impact melts. 
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Although X**® ( E In baDoltlo plagloelaoo waa 

racogniaed as a potentially useful crystallisation index {a,g. , 
Loaghi et al., 1^72| Crawford, 1973), the detailed roachanlama 
of Pe and Mg substitution in plagioclasa are debatable. Host -• 
chemical studies of lunar plagioclase indicate that nearly all 
deviations from Ab-Or-an chemistry can be explained by the 
components, Ca<Fe,Mg)SijOg and QSl^Og (e.g., Helll et al - , 1970| 
Monk and Hilda, 1973) where D signifies a Ca-site vacancy. 
Spectral studies (Bafaer et al, , 1971 and Appleman et al. , 1972} 
indicate that much of the Fa la not situated In tetrahedral 
sites as would be expected from a Ca{Fe,Mg)SijOg component, but 
appears to be situated In the larger, six- to eight-coordinated 
Ca-sites. Furthermore, Hafnor aj[. (1971} suggested that the 
Pe apparently in tetrahedral sites may be associated with 
oxygen vacancies. In order to predict the corapoBitlonal 
dependence of Fe and Mg substitution it is necessary to know 
the dominant Po-Mg component and also to' model tha-effects- of r 
substitution on two types of sites, 

ionghi et al . (1976} , via projections in the system 
Ca4SljOa-AlgSi_20g-si^Og-(Fo,Hg)2SijOg, showed that all the 
chemical variation in plagioclase from mate basal'ts can be 
explained by mixing the conventional plagleclaae components 
(An, AbandOr) with OSi408< 'Ca(Pe,Mg)Ei30j and (Fe,Hg)Al28i20g. 


H 

The data rule out any slanifleant substitution of oosponefits 
involving oxygen vacancies (e.g. , CaPnigSlgO^) , and also show 
that no more than about 10J( of the Fe and Mg in plagioclase 
from mare basalts substitutes In the form of the {Pe,Mg)Al2Sl20g 

component. Models of the multi-site crystal/single-site 
liquid equilibria (longhl, 1976) based upon the analysis of 
Orover and Orville (1969) show that deviations from one-site/ 
one-site equilibria will be obscured when the proportion of 
distinct sites in the multi-site phase is not 1:1. Thus In 
calcic plagioclase, where the proportion of sites available 
to Pe and Mg la likely to be l(ca}:2(Al), P« and Mg substitutions 
may effectively behave as if plagioclase had only a single typq 
of available site. 

We here examine the substitution of Fa and Mg in 
plagioclase. as a function of composition and temperature, 
recognising that the dominant component la Ca(Fe,Hg)SijOg. 
Possible pressure and kinetic effects will be considered. 
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Dlatrlbutlon Coefflelenta — Theory 

In order to deal with the oompoaltlonal dspcndenoe of 

the distribution coefflelenta, we will aaaume that all the 
+2 

Fe plus Mg forma, a CaCFe.MeJSl^Og component In plasloolaae. 
Consider the reaction: 

CaOj^ + PmOj^ + 3 aiOgj^ » CaFraSl^Ogp (j) 

with equilibrium constant: 


„Fm - ECaFmSl_On]„ 

eq — - I - ■ , # 2 % 

ECaO]j^ • CFm03j^ . CS102'3£ 

where Fm denotes either Fe or Mg, subscrlpta P and t. denote 
plagloclase and liquid respectively, and the brackets Indicate 
activities.- The conventional molar distribution coefficient Is: 

p,. (FraOJn _ 

* tSlO]^ • 

/ (3) 

where 0.2 Is the approximate stoichiometric factor for calcic 

plagloclase that relates PmO concentration to CaPmSl,0o 

3 8 

concentration when analyses are recalculated In terms of 
conventional oxides (l.e., AlgO^) and y Is an activity 
coefficient. Even If all the solutions were ideal mixtures of 
oxides the distribution coefficient would still vary as the 
concentration product, (CaO)j^* (SIO^)^. On the other hand. 

•If moat Pm substituted for Ca to form a FmAlgSlpOg component 
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then would vary Inversely with (Al 202 >j^» (SlOg)^. He may 
predict that with decreasing Ca/Al In the melt, the proportion 
of Pm substitution for Ca will increase. 

The direct composition dependence of simple distribution 
coefficients may be circumvented by the use of exchange 
distribution coefficients" such as; 

..Pe-Hg . <P«0)p* Kff,. 

■ rugoyj,* {PeO)£ “ 


,(Fe PeO CaMgSl Oft 
eg . 'h ^ Tp 30 
uMg ’ ,„MgO * “daPeSi-Oft 


Here the division cancels out the activity terms In equation 
(3). It can be readily shown that: 


.Pe-Mg‘ ^ 
‘P-L 


(i-xP«).x[« . 


(5) 


where X*'® = Pe/(Pe + Mg). Thus It la possible to relate X^‘ 
Pe ° 

to Xj^ by means of a single coefficient over a wide range 


of liquid composition. Given the formulation of plagloclase 


activities In terms of phase components, the activity 
coefficients and y^®'^®®^ 3°8 will probably vary with 

the An/flb ratio. Because of the division In equation (<l), the 
exchange coefficient, will be less sensitive to An/Ab 

variation than the single element coefficients. 


An alternative approach to avoiding composition dependence 
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la to define raope complex, sinsle-elemeafc distribution 
coefficients that explicitly wary with liquid compoaltlon 
(Drake and Weill, 1975) BUoh as 


. CCaO]j^. CPeOlj^- ( [SlOglj^)^ 

Several options are available for' approxlmatlns liquid 
component activities by assumlns Ideal solutions: 1) Ideal 

oxide mixing (activity equals concentration) 5 2) mixine of- 

Ions or' oxides on two separate lattices (Temkln, 19^*5)} or 
3) separata mixing among network-forming and network-modifying- 
components (Bottlnga and Weill, 1972). Although component 
activity approximations based upon the Bottlnga-Wetll melt 
model seem to- have* the greatest potential utility, these 
approximations can be used only with liquids with restricted 
ranges of composition because the rigid, norm-like calculation 
of the proportions of melt species leads to unreal disconti- 
nuities In the distribution coefficients. For example, in , 

an alkali-free liquid, if (CaO), - (A1„0„), then [CaO], - 0 

Po * * - 

^p-D ” “■ Therefore, we will adopt the ferakln -meit model 

as a tentative basis for activity approximations and calculate 

activities as follows: 

large cations (Fe, Kg, Ca, Cr, Tt, Ha, K) — 

[FeOjj^ ” ^Pe + Hg ,+ Ca + Ha T'K + ti + Cr^t 

etc. , and tetrahedrally coordinated Ions (SI, Al)' ~ 




^Si + Al^h 


where * signifies the Temkln model. 



All analyses ware made witli an AR£> SKK-SM 120000 electron 
microprobe at.lSkv, O.OIX to, 0.015' wA sample current and a 
1-2 pm wide beam. Synthetic, oryatalline ahorthite was the 
major element standard forCa, Si and Al and as the background 
standard for Ha, K; Fe, Mg and Tl. Typically, approximately 
1,000 counts were accumulated for Pe and Mg, giving a counting 
error of 10.03 mole STPeO and HgO. Relative counting errbra 
for Ca, Si and Al were < ±iS. ' 

Whenever the point of an analysis on plagloclase was, 
closer than 50' pm from a boundary with an Iron-bearing phase 
a correction was made for secondary flourescence of iron 
according to the empirical' relations 

®An-'®92 “ ^An''*w “ (-Q.037X),. 

determined by Longhl (1976) where C^® la the apparent concen- 
tration of Fe In anorthlte, I is x-ray Intensity, 92 refers 
to the Iron-rtch phase at x pm from the point of analysis. 

Due to the cubic dependence of secondary flourescence upon 
atomic number (Smith, 1965), the effect for magnesium Is 
negligible. Most eorreetians were in the range 0.05 to 0.20 
mole % PeO. Ho analysis vias performed upon crystals with 
minimum dlmenelona cio pm. 

Fourteen glasses with compositl.ons similar to .those of 
lunar basaltic liquids saturated with plagloclase were 
synthesized from reagent grade oxides and Amelia albite by 
multiple fusions In platinum crucibles under a PSXNg / 5S 
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atmosphere. After each fusion the sample was ground to a 
fine powder In a tungsten-oarblda ball mill. The composltiona 
spanned a rar^e of atomic Pe/{Fe + Mg) from 0.2 Co 0,8 and 
were chosen so that plagloclase was the sole llquldus, phaao,. 
the llquldus temperature was approximately 1200“0, and the 
concentration of TlOg was approximately 1 wtj!. 

Experimental charges wore loaded Into hlgh-purlty iron 
capsules (Walker et al, , 1976) enclosed In sealed, evacuated 
silica tubes, fhe silica tubes were suspended in a Oeltec 
T31 furnace with MoSlg elements and dropped into water at the 
end of the experiments. The charges were held ''<20°C above 
the liquldus for a mlnimum.of 4B hours to ensure equilibration, 
with the iron capsules, then the temperature was lowered five 
bo twenty degrees below the llquidus in one-half hour and 
held at constant temperature for a mlnlraun of ona.-week. , This 
technique allows plagioclase crystals to grow with minimum 

dimensions > 10 un and to approach equilibrium with the melt 

« 

(Drake and Weill, 1975). 

Experimental products of the seven compositions containing 
alkalies were plagued by the presence of Iron spherules inoluded 
in the plagioclase. Because of the uncertainties Involved in 
applying corrections for secondary f lourescence , the data for 
these experiments (honghl, 1976) are not Inoluded in the present 
discussion, although the distribution coefficients for Iron and 
fliagnesiunt are aimllor to those obtained from tho alkali-free 
compositions. 
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A series of five runs of variable length were made at; 
1230®C with a single composition, DPC.8, in order to test for 
possible Iron gain as a result of oxygon diffusion through 
the silica glass tubing (Kesson, 1975). There was approxi- 
mately 0.5 wt« PeO gain in the first t8 hours, but no 
measurable variation • from liB to 2t0 hours. Thus oxygen 
diffusion la not a factor In these experiments. However, 
changes in the composition of the melt as it equilibrates 
with tho capsule are a potential source of serious error in 
short term experiments. 



11 



The distribution of Iron and magnesium between plagloalase 
and liquid can be measured directly fcy mlcroproba analysis' of 
quenched natural and synthetic oryotaUlquld pairs. Indirect 
Information from crystalline rocks can be obtained by comparing 
Fe and Mg concentrations In the cores of the eBrlleBt-forroed 
plagloclaae crystals with the estimated composition of coexisting 
liquid. If plagioclase was a llquldus phase of a quickly cooled 
melt, than the bulk rock composition represents the coexisting 
liquid composition (e.g., 68lll5, ItHO, 12038, 75035). If 
plagloelase crystallized late, then the coexietlng liquid 
composition may be estimated from malting experiments If the 
rook has cooled slowly (e.g., 15065, 15555 and 7001?). Finally, 
if the rock has homogeneous -minerals (e.g., pink-spinel and 
trootollte) then the FeO and MgO concentrations in. the. coexisting 
liquid can be calculated from olivine compositions if approximate 
temperatures are known (e.g., Boeder and Emslie, 1970). 

Tables 1 and 2 contain direct data on eleven sets of 
synthetic and eight sets of natural plagioclaoe liquid pairs. 

The natural samples include aubraarine baaalt glasses and 
glasses from Makaopuhi lava lake (Hawaii). Estimates of liquid 
composition are made for eight lunar feldspathic and mare 
basalts, and for three terrestrial and one lunar* plutonlc rock • 
(Tables 3 and <t). The various distribution coefficients are 
listed In Table 5, 
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Baaaltio gamnlea , 

Figure la shows a close approach to .04 for both 

lunar and terrestrial baaaltio samples. Recorded MgO ranges 
up to 0.75 mole S, The hlgh-Tl lunar basalts show the great 
deviations. However, Figure lb shows that terrestrial 
plagloolase has significantly higher iron concentration than 
lunar pXagloclase at a given concentration of iron in the 


liquid. Given the similar- values of terrestrial and lunar 


the major cause of the higher Fe concentrations in terrestrial 
plagloolase is not likely to be an effect of the lower 



ratios of the terrestrial plagloclase, but rather 


the presence of ferric iron In terrestrial plagloclase-iiquld 


pairs. 

t'Po" 

"P-L (where "Pe” » Fe"*^^ + Pe^r) to be greater than 

(Fe'*’VFe'*^®)p must exceed (Fe*^/Pe’*^^)j^. Assuming that 
Pe*^^ 

values of Kp_^ ore similar for lunar and terrestrial crystal- 
liquid pairs as are the values for between two-thirds 

and three-fourths of the Iron in the terreatrlai plagioclase 


is ferric. 


Figure 2 shows the exchange relations of Pe and Mg for ' 
the natural plagloclase-llquld pairs. The value of 

r*-L» 

('u’o.5) for the lunar low-Tl basalts means that at low oxygen • 

fugaoity Xp® will be considerably lower than X^®. 

increases with Increasing oxygen fugaclty and has an average 
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*‘Fe** J'Fe*' 

value of l.ft foi> terjvsfltrlal basalts, ao Xp > Xj^ . 
Presumably, if we could correct for all the ferric iron in the 
analyses of the terrestrial plagioclase and glaasea we would 
observe an average value of close to 0.5. The. 

variation of with oxygen fugaolfcy provides a crude 

oxygen 'barometers In environments nioi*e oxidising than those 
of terrestrial basalts, valuea higher than 

l.t. 

Synthetic* Plagloelase^Moui'd Pairs 

Figure 3a shows the simple Hg distribution between . 
synthetic plagioclase and coexisting liquid. Even though the 
temperatures of equilibration spanned a limited range <1218*C" 
1189®C} as opposed to the range of to 1050®C for the 

natural samples, there la a good deal less oohorcnca to the 
distribution between synthetic crystal-liquid pairs! the 
relative error (= ±l 0 /Kp_j^) for the terrestrial Kp®j^ is ±22j{, 
whereas for the synthetic pairs it is ±30J5. In Plgure 3b we 
have plotted the HgO concentration in plagioclase versus the 
Temkln activity product [CaO]^* [MgOl»* (CSlOglj^)^. The 
relative error of the average is now *15?. Thus much of 

the variability of Kp®j^ in Figure 3a is due to compositional 
effects. The plagloolase-liquld distribution of Fe shows a 
similar compositional dependence although the reduction of the 
relative error is not as greats “ 0.033 * 0.010 (*30?), 

whereas » 0.0?2 t. 0.017 (*2hJO. Errors from Inaccurate 

application of the Fa secondary flourescence correction at 


It 

Incllneil plagloolase-glasa boundaries may cause higher relative 
errors in than Kpfj^. 

Besides -lack of complete equilibration and multi-site 
substitution, another factor which may affect • the- Inferred 
magnitude of the single element distribution coefficients is the 
diffusion controlled buildup' of rejected elements marginal to 
grotjlng crystals (Albarede and Bottinga, 1972). Because the 
plagioolase-liquld distribution ooefffeienta of Ps' and Mg are- 
much lesa than one, they are particularly susceptible to errors 
. Induced by the under-estlroatlon of FeO and MgO contents In the 
liquid. Longhl (197fi) reported negligible enrichments of HgO 
Bsarginal to plagioclase from Makaopuhl lava lake, ao relatively 
slowly cooled naltural plagioclase apparently does not yield 
erroneously high distribution coefficients. However, Kirkpatrick 
et ar. (1976) report more than two orders of magnitude differ- 
ence in the growth rate of anorthlte at 5* to 20®C undercooling. 
Thus it la reasonable to expect some variability in the inagnt- 
tudeo of the single element distribution coefficients among 
oynthotlc crystal-liquid pairs grown at different undercoolings, 
and ainoo these undercoolings are probably larger than those of 
the more slowly cooled natural liquids we can expect synthetic, 
single element distribution coefficients to be somewhat larger 
than those determined from the more slowly cooled natural ^ 

samples (Table 5). Since, however, ®hd Kpfj; have similar 
magnitudes, there will be little kinetic fractionation of Pe 
and Mg, and thus natural and synthetic values of ought 

to be similar {Table >51. 
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Figure 4 shous the Fe-% exchange relatione for the 
synthetic plagloclaae-llquld pairs. Much of the spread In 
data points is prohably duo to errors in applying the eorreotlon 
for Fe secondary flourescenoe. The 'relative counting errors do 
not exceed ilOjC of The average value of 

of the natural low-Tl basalts ,<• 4 9) is. well within the error 
range of the synthetic Kpf£”® (.51 * .13). However, the 
average for hlgh-Tl basalts of ,33 is slgnlfteantly 

lower. ' Longhl- et al. (1975) noted a similar effect of T1 upon 
^01-^ ■ shown In Figure 4 Is the curve of constant 

*^01-L* *’ ’determined by Longhl et (1975) for lunar 
basaltic liquids. We see by comparison that at low oxygen 
fugacity and low total pressure the iron to magnesium ratio 
of plagloclase will be lower than that of coexisting low-Tt 
liquid but higher than that of coexisting olivine and 
probably orthopyroxene. Plagloclase ooexistirig with a highrTi 
liquid will have similar to coexisting olivine and 
orthopyroxene. 

The absence of any lar^, systematic deviations of the 
data' points in Figure 4 from the average value of 
Indicates that effects due to the two-site substitution of Fe 
and Mg are small in the range of from .2 to .7. 

Plutonic Samples 

The limited analysis, of plutonlc plagloclase-llqutd pairs 
Is based upon three terrestrial samples (two from the Skaergaard 

“ I 

intrusion — lower zone, one from the Ardnanurohan eucrlte) and 


one lunar — pink spinel troetollte (Prlnz et a^. , 1973). Liquid 
concentrations of PeO and %0 are calculated from the compositions 
of coexisting olivines by the equations of Boeder and Ensile 
(1972) and Longhl et (1975) by estimating temperatures of 
equilibration. For the terreatrlal liquids PeO is converted 
to "PeO" (« PeO + PeOj j.) by assuming (Pe^VPe'*^^)j^ ■« 0,2 
(Longhl, 1976), Plutonic plagloclase-llquid pairs have 
considerably lower values of single element distribution 
coefficients, but higher values of Kp*£ than do their low 
pressure counterparts (Table 5). The first observation is the 
result of lower ooncentratlons of FeO and MgO in plutonlc 
plagloclase than in basaltic plagleclase. As noted above, the 
absence of significant HgO gradients in the melt marginal to , 
lava lake plagloclase precludes a kinetlc/diffuslon explanation 
for the difference In Xp®j^ and Kp®j_ between lava-lake basalts 
and plutonlc rocks. In order for the differences in distribu- 
tion, coefficients to.be due entirely to a pressure effect, the 
partial molar volumes of CaPeSi^Og and CaMgSl^Og must not only 
be considerably larger than that of anorthlte, but also be 
approximately equal to one another. The latter is doubtful. 

An alternative explanation of the differences in 
distribution oocfflclonts lies in the dependence of the 
distribution coefficients upon liquid undercooling (or super- 
saturation). Hopper and Uhlmann (1974) show that the kinetic ’ 
crystal-liquid distribution coefficient _ for a growing crystal 
la a function of undercooling, and that for Isothermal growth 
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tue compositions of crystal and liquid are closer together than 
the equilibrium compos Itlona. Lofgren (1973) has demonstrated 
this effect for plagioclase. The smaller the undercooling . 
the closer the approach to the equilibrium compositions, thus 
plagioclase growing under oondltiona of alow cooling and hence 
relatively small undercooling’ le able to exclude lron> and 
magnesium more effectively. More rapidly cooled basaltic 
plagioclase shows- the kinetic* effects of larger undercooling 
with a consequently enhanced Incorporation of Iron and 
magnesium. 

Whether this kinetic effect or the pressure effect can. 
adequately explain the differences in observed distribution 
coefficients among plutonlc and volcanic rooks is not known. 

It seems best, however, to treat distribution coefficients 
from these two types of rocks separately at this time. 
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Thermal Dependenoe of the Distribution Coefflalenta 


The similarity of values of for natural plagloclaae- 
liquid pairs that span a range of temperature of more than 
300‘’C is an Indication that the temperature dependence Is 
small, the temperatures of the submarine basalt liquids have 
been calculated from ollvlne-llquld equilibria by Prey et al_. 
(197^); the temperatures of the lava lake samples are obtained 
by direct measur^ant and are Hated by Kirkpatrick (1976). 

Log Kp®j^ and log for the terrestrial basaltic samples are 

assumed to be linear functions of 1/T“K. Regression equations 
In the form log K » fl/T®K + B are listed in Table 6 along with 
similar equations for olivine determined by Longhl ^ al. (1975), 
The. ollvlne-llquld dlatributlon coefficients, and 

show similar temperature dependences; It la likely that the 

jpA 'PaH 

temperature dependencies of Kp_j^ and will be similar to 
Kpfj^ and Kp*^ respectively. The plagtoolase-llquld pairs have 
smaller values of A Indicating a -weaker thermal dependence than 
the olivine equilibria; also the linear plots are relatively 
poor fits to the plagloolase data. Therefore, neither log 
nor, log Kp®J^ is a useful geothermometer; on the contrary, single 
values of the plagloclaoe-llquid distribution coefficients may 
he employed throughout the crystallization interval with little 
thermal error. We can expect the temperature dependence of 
to be even smaller than that of Kpfj^. 
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Applteatlona 

Since Pe and Mg aubstltuta in a aystematio, preilctablo 
wanner In plagloolaae we aay' employ tha ooncentpatlono of 
these elements in plagloclase to calculate liquid compositions 
and as a oryatalUzatlon Index, A crude, but potentially 
useful application to oxygen bayometry-haa been noted above; 

And now that distribution coefficients are known it may become 
possible to calculate from tha models of Smith et al. (1955) 
end Hopper -and Uhlaann (197t) the crystal growth rate of 
plagloclase In quenched melts by measuring the diffusion 
controlled gradient of MgO in the melt marginai to plagloclase. 
There are three simple applications to lunar petrology. 

4 

Anorthosites 

The presence of plagloclase In tha anorthosites, with 
relatively high anorthlte contents (An’gg^gg) coupled with 
pyroxenes of intermediate (0.3-0. 6) poses a thorny 
interpretive problem. Smith and Steele (197^ )» noting that 
much of the pyroxene in these rooks In the form of rods and 
stringers a few microns wide appears to be exsolved from the 
plagloclase, suggested that of these pyroxene crystals 
represented the original X**® of the plagloclase and that, 
therefore, the of the liquid from which the anorthosites 
crystallized could be calculated from They 

obtained a value' of = 2 from three published 

plagloclase-llquld pairs and were thus able to calculate 


ffa 

values of XJ” leas than O.a, uhieh in turn led them to the 
conclusion that the anorthosites had crystallized from a 
relatively primitive liquid. In view of the dlBtx^lbutlon data 
presented above ouch an origin la unlikely. 

Since no corrections for Fe secondary flourescence were 
made in the original work, we are certain that the value of 

“ 2 Is erroneously high, Olven the measured values'Of’ 

c 1 for lunar basaltic and plutonlc rocks, we may con- 
clude that if the small rods and stringers of pyroxene are in 
fact the result of the complete exsolution of Pe and Kg In 
plagloclase, then the anorthosites crystallized liquids with 
intermediate X^®. 

The application of distribution data to anorthosites is 
by no means simple, however, because Walker et al^. (1973) and 
Dixon and Fapike (1975) have shown that plagloclase in the 
anorthosites (e.g., 60025) may retain measurable amounts of 
Pe and Mg in the presence of these small pyroxenes. If these 
small pyroxenes are exsolutlons then neither the X^® of 
pyroxenes nor plagloclase may represent the original Xp®. If 
tha small pyroxenes are primary inclusions then some metamorphlc 
exchange of Fe and Kg between pla’gloolase and pyroxene may 
have altered the original x^® of both pyroxene and plagloclase. 
Regardless of wblch of these possibilities la correct the 

orleinal Xp In 60025 muat be greeter than or equal to the 
Fe 

present a 0,35 0,t5 and less than or equal to the present 



21 


X^® > 0.5 (Walker et al., 1973). Application of a lunar 
(0.5 - 0.7) then leads us to the oonclusloii that 60025 
crystallized from a liquid with Pe/(Pe + Mg) “ 0.6 - 0.7, 


Zonatlon 

There are many published reports of Pe and Mg zonatlon in 
basaltic plagloclase, however, some workers have been skeptical 
of the conclusions reached because the ■ secondary flourescence 
effect alone will produce an apparent Pe zonatlon. We, 
therefore, present several patterns of flourescence-correoted 
plagloclase zonatlon to demonstrate the utility of Xp 
variation. 

In Figure 5a, "anorthlte'* content ( gg ' ^a T ~K^ plotted 
Fe 

as a function of Xp for plagloclase cores, mantles and rims In 
two slowly cooled mare basalts, 15555 and 15065, with plaglo- 
clase poiJcllltlc texture. As Xp® increases, anorthite content 
decreases much as expected. Zoning in most Individual crystals 
follows the general pattern— cores have the highest ga t ' Wa V X 

the lowest X^®, whereas the rims hav« the lowest ca ‘ " ' t " ' Sa '~ + ’ " lt 
and highest X^®— an indication of outward growth. However, note 
that one of the crystals In 15065 has Its highest ca' t Ha T K 
and lowest X^® at Its rims. This apparent reversed zoning is 
consistent with nuclcatlon of the plagloclase upon separate 
pyroxene crystals and subsequent growth Inward Into the volume 
between the pyroxene crystals. 

In Figure 5b are similar plots of plagloclase conation for 
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two-flnar grained msre basalts, 12002 and 12021, in which 

plagloclase laths oomzwnly have hollow cores ^ In contrast 

to Figure Sa there Is a small, overall increase In anorthite 

Fa 

content with Increasing X In the plagloclase. cral^ford 
(1973) has previously observed this phenomenon in 12021, Tblo 
sort of reversed zoning Is consistent with delayed nucleatlon 
of plagloclase, followed by plagloclase crystallization at 
progressively decreasing under-cooling at near Isothermal 
conditions (Hopper and Ohlmann, 197A{ lofgren, 1973). 


» 


Variation of Iron and roagnealuih In plagloclase may help 


Identify chemical heteri%eneity In Igneous rocks, a clue to 
an Impact origin (Slmonda, 1975). Figure 6 la a plot of 

Ca + Na + K Pe^+ Mg cores of the larger, texturaily 

primitive blooky crystals and laths In rook HSIO, a feld- 
spathio basalt. The variation of t Ma 'TTt ®°hetant 

is nearly twice as great as that in the mare basalts 
(Figure 5).. Further, there la a structure to the variations— 
cores with ca 1- Ha t K concentrations of KgO > .10 

mole g (generally greater than .15 i with a maximum of .25 JS); 

Ca 

cores with, Ca"" j^ ~ Ma ~ +~ K ^ have concentrations ' of KgO < .10 
mole i (generally less than .05 %)-. Such variation of 
c ' a'i Ha t K » F e ^~ Hg ^2° observed In plagloclase 

cores in the mare basalts which are truly volcanic. Although 
this rock was completely or nearly completely molten at one time. 
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1*1310 crystalljzed fpom a nsolt that was inhomog'enaaus . (toxygan 
at al , ( 1972 ) Identifieii rock 1^310 as a probably impact melt 
on the basis of high, non-lunar oidoropblle element abundances, 
Schonfleld and Ksyer (1972) showed that the bulk chemistry of 
this rook can bo explained largely as a mixture of KfiEEP and 
anorthosltlo' components. It seems logical to conclude, there- 
fore, that the plagioclase cores with low j( and 

hlgh-K are-ghosts of KREEP materials incompletely homogeniaed 
during the melting process. 
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Figure la. Simple molar Mgo distribution between calcio 
plagioclasQ and natural basaltic liquids. Lunar liquid com- 
positions are all estimated; 

lb. simple molar "FeO" distribution between calcic 
plagioclase and natural basaltic liquids', where "Fo0“ is total 
iron as FeO. Symbols same as in la. 

Figure 2. Exchange relations of Fe and Mg between plagioclase 
and natural basaltic liquids. Symbols same as in Figure 1. 
Figure 3a. Simple molar. Mgo distribution between synthetic 
plagloclase-llquid pairs. Temperature range is 1189» to 121B“C. 
P = pyroxene present. 

3b. MgO concentration in synthetic plagioclase as a 
function of Tembin model activity product; 

Figure 4 . Exchange relations of Fe and Mg .^between synthatlor ■ 
plagioclase and liquid. from Longhi et al . (1975). 

Figure 5a. Plagioclase zonatlon in coarse grained, mare, basalts,’ 
1555S and 1S06S with plagioclase poikilitlc texture. C ^ core, 

R “ rim, 

Sb. Plagioclase zonatlon in fine-grained mare basalts, 
12002 and 12021, with hollow plagioclase laths. 

Figure 6. Anorthite content versus (*^®/Fe + Hg)j, in plagioclase 
cores from lunar feldspathic basalt, 14310. High K ff mole % 

KjO > .10. 



'able 1. PlaRloclsse - liquid: distribution data: Terrestrial basalts 



Table 2. Pla^ioclafie-liguid distribution data 

t 

Synthetic Alkali-Free coraposltiona) Fe Capaulea 



RFC.2-18 

RFC. 3-18 

TFC.4-12 

TFC, 4-4 

Qrc.S-17 

DFC.6-13 

Wt4 SiOj 

54.18 

54.77 

57.76 

56,15 

52.78 

51.85 

TiOj 

.90 

.87 

1.07 

1,11 

n.d. 

1.16 

AI 2 O 3 

15.59 

14.81 

16.13 

15.05 

13.77 

15.19 


0 , 00 ' 

.11 

.14 

n.d. 

n.d'. 

. 22 ' 

FeO 

7.68 

7.68 

4.87 

6.19 

8.12 

12.74 

HgO 

8.32 

7.94 

9.24 

3,58 

6.73 

5.89 

CaO 

11.75 

12.27 

11.98 

11.75 

15.66 

12.65 

KjO 

0.00 

O.OO 

.01 

n.d. 

n.d. 

.01 

KajO 

0.00 

.07 

O.OQ 

n.d. 

n.d. 

0,00 

Sum 

98.42 

98.53 

101.13 

98.82 

97.06 

99.71 

atomic 

Fo + Mg 

.340 

.352 

.220 

.288 

.404 

.548 

Temp (”C) 

1224/1109 1224/1189 1257/1210 1241/1212 1233/1199 1257/1213 

■ Time (hra) 

48/168 

48/168 

SO/240 

48/210 

96/225 

50/240 

1 

Other phases 

opx? 

opx 

- 

opx 

> 


avg. mole% 







<PeO)p 

.26 

.26 

.12 

.24 • 

.32 

.36 

(MgO)p 

1.06 

1.03 

.84 

.85 

l.DO 

.60 

plag. analyzed 

2 

• '3 

2 

5 

S 

2 


Table 2 (continued). Pl& 3 ioclaBe-|i^uld diatrlbutlon data 
Synthetic Alkali-Free Composlilonsj Fe Capaules 


■ 

DFC.e-14 

OFC.6-15 

DFC.7-3: 

2 DFC.7-15 

CFC.8-18 

wtt SiOj 

49.94 

51.16 

45.79 

43.48 

49.13 

TiOj 

1.25 

1.25 

1.06 

1.12 

1.18 


14.69 

14.13 

14.15 

13.71 

13.16 


n.d. 

n.d. 

n.d. 

n.d. 

n.d. 

FeO 

12.17 

12.35 

19.20 

16,31 

16.05 

HgO, 

6.29 

6.50 

S.79 

6.2? 

S.47 

CaO 

12.59 

12. 4S 

11.69 

11.99 

12.30 

KjO 

R.da 

n.d. 

n.d. 

n.d. 

n.d. 

HajO 

n.d. 

n.d. 

n.d. 

n.d. 

n.d. 

Sum 

96;93 

97.03 

97.68 

97.88 

97.28 

atomic 

.520 

.516 

.650 

.593 

.622 

Temp (”0) 

1240/1217 

1241/1212 1241/12: 

.2 1233/1206 

1233/1206 

Time (hrs) 

40/211 

48/210 

48/2K 

II 48/168 

48/168 

Other phases 

- 

- 

- 

- 

- 

avg. mole! 






(FaO)p 

.39 

.34 

.44 

,35 

.31 

♦ (MgO)p 

.60 

.66 

.47 

.53 

.39 


no. plag. analyzed 4 $ 8 S 3 



TABLE 3 

PLAGIOGLASE-LIQUIC DtSTRIBUTIOS DATA: NATWHAL LUSAR BASALTS 


Wt»-L 


feldspathlc 


low-T± mare 


hlgh-Tl mare 

681115^ 

68415 - 
(102 An^Q)* 

s 14310^ 

15555^ 

15065^ 

12038^ 

70017^ 

75035"® 

SlOj 

«5.3 

45.7 

48,27 

45.8 

46.1 

46.56 

41.1 

41.6 

TiOg 

.35 

,39 

1.27 

3.26 

2.53 

3.31 

10.3 

10.2 

AX^O^ 

28.7 

28.0 

20.26 

10.8 

lX.l 

12.53 

9.99 

9.75 

crgOs 

.14 

.16 

.20 

.29 

.31 

.27 

r-i 

• 

,27 

PeO 

3.69 

4.12 

6.11 

20.0 

20.1 

17.99 

19.0 

l8.5 

HgO 

4.44 

4.96 

7.76 

5.74 

5.55 

6.71 

6,44 

6.16 

MnO 

- 


— 

.31 

.32 

.27 

- 

.30 

CaO 

16.3 

16.0 

12.25 

11.7 

ii.4 

11.62 

11.3 

11.5 

KgO 

.12 

,13 

.81 

- 

- 

.07 

.05 

- 

NSgO 

.5 

.54 

-55 

.14 

• .17 

.66 

.20 

.43 

Sum 

100.4 - 

100.0 

99.48 

98.0 

97.7 

99.99 

98.7 

98.8 

atomic 









Fe/Pe + «S 

.31 

.31 

• .34 ■ 

.66 

.67 

.59 

.62 

.63 

T°C 

.1,1425® 

<1400“C 

4,1325° 

1150° 

1147° 

1160°^ 

1144° 

1137° 

avg. mole % 









PeOp 

.04 

.08 

.16 

.40 

.35 

.25 

.32 • 

'.34 

MgOp 

.20 

.35 

,48 

.38 

.30 

.50 

.68 

.54 

•sample may 

contain laineited or xeno crystal 

material- 

-An^g 





X. Walker et al. U973J 5. Longhl et al. (X97'i) 

2. Walker et aX. (1972) 6. Walker et aX. (X975) 

3. unpublliHe3“experiBentaX 6ata 7, Blgsar et IT. (19713 

A. Conpstcm ^ (1971} 


TABIS 4 

PLAGIOCLASE-LIQUID DISTRIBUTION DATA: PLUTONIC ROCKS 


lunar terrestrial 


Hole S 

Pinkr-spinel troctolite^ 
6?ti35 65785 

Skaergaard 
. L2A 

Skaergaard 

LS.-ffflC 

Ardnamurchan 

eucrite 


5.0 

10.8 

24.7 

17.2 

19.6 

»8°01 

60.9 

55.8 

41.9. 

46.7 

■46.6 

PeOj,2 

5.0 

10.8 

17.3 

12.0 

13.7 


20.0 

18.4 

6,91 

10.4 

9.99 

npgci«^ 

- 

- 

20.8 

14.4 

16.4 

T°C 

1270°^ 

1200°** 

i( 

1200° 

1200®** 

"PeO”p 


.05 

,33 

.29 

.63 

MsOp 


.17 

.06 

.04 

.08 


X. Olivine preoeael plagloclase in 6?il35> tut crystalliaea after plagloolase in 65785* 
heaoe the actual KeOj^ and are taken to he the average of the computed values. 

2. Calculations based upon equations of Boeder ai^ Emslle (1970) and Longhi et al. (1975)' 
see table 6 of this study. 

3. Walker etal. (1973). 

4. Estimated. 
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CISTRIBUTIOH COSa^FJCIENTS 


Wole Praotilon 



„.Fa 

*P-t 

j,Mg' 

“P-L 

„Pe» 

“PriL 

j,Hg» 

‘'P-L 

vPe-Kg 

^P-L 

Lunar 

. 





Basaltic 






loM-Tl 

.030 

.041 

.056 

.12 

.49 

hlgh-Tl 

.020 

.061 

.075 

.24 

,33 

Plutonic 






(P.S.-P.) 

.006 

.009 

.03 

.04 

.7 

Terrea trial 

1 





Basaltic 

.osa 

.041 

.18 

.13 

1.4 


+.008» 

+.009 

+.oa 

+ .03 

+•5 

Plutonic 

.025 

.0063 



4,2 

Synthetic 

.033 

.065 

.072 

.14 

.51 


+ .010 

+.016 

+ .017 

+.02 

+.13 


, TABtE 6 

'tEMPERATlffiE DEPENDENCE OP DISTRIBH'PION COEPPICIKHTS 

log K »• + Bj r »• correlation coefficient 


K 

A 

B 

r 

rRS 

'^P-L 

1216 

2.265 

.48 

^P-L 

1255 

1.787 

.50 

„Pe * 

2818 

1,212 

.97 

„M8 ‘ 

2616 

1.625 

.35 


'longhl al.. (1975)— lOH-Tl ollvino-llquid palra 
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PLAGIOCLASE-LIQUID PAIRS 



MOLE% M60-Lia 


» I 

P/nji 








SYNTHETIC 



SYNTHETIC 





Lomkt ^laX. 

• K'buhf ■ 


TS^s pm ® 5>i)ir/; 



fhf H 
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ABSTRACT 

The Apollo 12 olivine basalt suite shows a strong positive 
correlation of grain' size with normative olivine content. 

This correlation is interpreted to mean that the suite of samples 
represents the basal portion of a cooling unit which differentiated 
by simple olivine settling. The grain size of plagloolase 
observed in the coarsest samples indicates the cooling unit 
may have, been as much as 30 meters thick. The amount of 
olivine concentration observed in the suite is quantitatively 
Internally consistent with simple olivine settling in a 
magma body of this size which has the composition of the chill 
margin. 

INTRODUCTION 

The Apollo 12 olivine basalts form a coherent group whose 
major element chemistry can be related by olivine fractionation. 
Sufficient petrography and experimental petrology Information 
exists to make a detailed analysis of the differentiation process 
of the magma body (or bodies) in which this suite of samples 
formed. This analysis is based on a simple model of gravltatlve 
olivine settling. Quantitative input for the model calculations 
comes from crystallization experiments at controlled cooling 
rates on one of the members of the suite . ' These 
experiments quantify the cooling history and the physical 
properties of the magma during solidification. It is the 
purpose of this paper to show the relation between the cooling 
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history, aifferentlation, and petrographic features of this 
* 

suite of olivine basalts. Normative olivine abundanoe Is the 
principle index of differentiation and grain size Is the prin- 
ciple feature related to cooling history in our dltousslon, 

30 our task is to explain the strong positive correlation 
between olivine abundance and grain size which is observed 
in this suite. 

CHEMISTRY AND RETHOORA?HY 

'A sub-group of the Apollo 12 samples~tbe oll\lne basalts, 
plcrites, and gabbros — has a major-elemenb chemical variation 
which can be described to a good approximation by an addition or 
subtraction of olivine of about (Kuahlro and Haramura 1971} 
Conpston et al. , 1971; James and Wright, 1972). Green et al. 

{1971) showed that the liquidus olivine of the least olivine- 
normative member of the group (l.e,, 12009) was of about this 
composition whereas the liquidus olivine of the more olivine- 
normative members was significantly more magnesian than any olivines 
actually observed in the suite. This led to the oo'icluslon that 
the less olivine-normative members of the suite repjresented 
liquid magma compositions into which liquidus olivine accumulated 
toform the rest of the suite — a conclusion compatlb|le with the 
observation that 12009 (the proposed magma composition) had a 
vltrophyrlo texture indicating quenching from a melji. Drove 
et al, (1973) and Walker et al. (1975) also supported this 
conclusion from work on 12002, a pierite near In cojsposltion 
to 12009. 
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Table 1 lists the normative olivine abundance for rocks 
in this suite taken from published chemical analyses. Hormatlve 
olivine varies by more than a factor of two. Modal data collected 
on 1 to 3 sections of each rock are also listed. The correspond- 
ence of modal and normative olivine is acceptable for samples 
12009', 12075, 12002, 12040, and 12035 but Is quite poor for 12009, 
12020, and 12018 which show too little modal olivine relative to 
the norm, msldentifloatlon of pyroxene for olivine is not 
thought to be the cause of the discrepancy because James and 
Wright (1972, Table 13) also record this discrepancy for 12018 
but not for 12020, This suggests that there may be significant 
heterogeneities in the distribution of the olivine phenoorysts 
in these samples. This is also suggested by the large range 
in the chemical mode for 12009 in Table 2 A of James and Wright. 

Figure 1 shows a series of photomicrographs of textures of 
rocks in this suite. The photos ave arranged in order of 
increasing plagloclaae grain sise and a very obvious relation 
to normative olivine content is apparent. The coarser samples 
are more olivine-normative. Sample coarseness IS determined by 
the grain size of the groundmass minerals. Although the olivine 
phenoorysts show a very limited range of size, there is a system- 
atic variation in their habit. In 12009, the vitrophyre, the 
olivine phenoorysts are euhedral to euhedral-skeletal. As the 
samples coarsen the olivines become subhedral and lose all 
distinction as phenoorysts in the coarse-grained, hoTocrystalllne 
gabbros . Pyroxene phenoorysts become equally Indistinct in the 
coarse-grained gabbros and they are also absent in the vitrophyre. 
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Since pyroxenes are present both as phenoeryst and grountoass 
minerals and therefore show quite a variety of Hiorpholcgles 
and sizes in a single sample, groundmass plagloclase is pre- 
ferred as an index of sample coarseness. Table 1 lists the 
average width of the largest plagloclase laths observed In thin 
section. The size varies by almost three orders of magnitude 
and it can be seen that there Is an excellent correlation of 
plagloclase slae with normative olivine, content. This is shown 
graphically in Figure 2. Plagloclase size is a useful parameter 
of coarseness since it is more easily quantified than pyroxene 
size. [Furthermore, groundmass minerals grown after magma 
emplacement, such as plagloclase, are not as likely (as Intra- 
tellurlo phenocrysts are) to be distributed according to size 
during magma flowj and hence they are probably a better Index of 
cooling history after emplacement.! The two minerals give gener- 
ally similar ordering of coarseness by "eyeball" as can* be seen 
from Che study by Hodges and Paplke (1975) who did nob use 
plagloolaae size specifically. Their coarseness ordering of 5 
samples reversed only the adjacent pair 13075-1200^1 relative 
to the present sequence of 3 samples. 

The fact that the chemistry of the samples can be related 
by olivine addition is assumed to mean that the samples are 
genetically related in this way. The specific composition of the 
magma will be taken up below. The attempt to reconstruct this 
suite of samples within a single magma body runs counter to some 
current thinking about -lunar petrogenesls . Previous petrographic 
studies by James and Wright (1972) and Warner (1971) have tended 


to separate the coarser and finer grained samples with the Impli- 
cation that they have a separate genesis. An exception to t^ls 
approach is the study by Hodges and Paplke (1975). It is per- 
haps the evidence from minor elements which in the past has most 
strongly suggested that any simple story implied by the correl- 
ation of petrography and major element chemistry is misleading. 

For example, Coropston et. al . (1971) showed that olivine fraction- 
ation alone cannot relate all the samples since at least. three 
different subgroups of the suite can be recognized on the basis 
of Rb/3r arid initial ®^Sr/®^Sr. Furthermore, H. A. Schmitt 
(personal communication) has pointed out that at least three 
subgroups, which cannot be related by simple olivine fractiona- 
tion, are recognizable on the basis of ratio. These group- 

ings are ndt the same as the Hb/Sr ones. It would appear that 
these arguments are very powerful and that It may be necessary 
to abandon the notion that the whole suite is part of the same 
cooling unit. However, it should be noted that substantial 
isotopic Inhomogeneity has been reported in Sr in a single recent 
terrestrial basalt flow (Laughlln et al. , 1972) which is pre- 
sumably related to cryptic contamination. If such processes 
supplement olivine fractionation as a control on minor-element 
and isotopic chemistry in the lunar ease, then it may net be 
necessary to abandon the idea that the samples are consangulnoua . 
It is worth noting that even if the samples are not consanguinous , 
the process we shall explore here would have to occur in each of 
however many magma bodies of the same major-element composition 


are required by the minor-element composition. Separate bodies 
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of different major-element composition would not provide a 
plausible explanation of the correlation of petrography and chem- 
istry since the simple gravltative settling process Is quite 
sensitive to magma composition. Although we cannot exclude the 
possibility that some other process may be responsible for the 
correlation we intend to ahow that simple gravltative settling 
provides a quantitatively Internally consistent explanation. 

CRYSTAL SETTLING 

We have done an experimental study of the crystallization of 
composition 12002 at controlled cooling rates (Walker et al,, 
1976 a). It was possible to specify the cooling history of sample 
12002 at various points In its crystallization by comparison of 
petrographic and experimental olivine nucleatlon density, pyroxene 
"swallowtail" spacing, and groundmass grain size. We concluded 
that 12002 cooled through its llquldus (-1325®C) at -0.8°C/hr. 
and that the cooling rate subsequently decreased, so that about 
2 months elapsed before the rock reached 1000°C. Prom this exper- 
mental work, the composition of the residual liquid and the degree 
of crystallinity are known as a function of temperature as are the 
sizes, shapes, and compositions of crystals which grow. Prom this 
Information, viscosities, densities, and crystal settling veloci- 
ties may be calculated as functions of temperature. This calcula- 
tion Is more practicable In the lunar case than In common terres- 
trial circumstances since there is no HgO present and negligible 
PGgOg. The presence of either can have a considerable effect on 
the result, and their values can only be estimated In terrestrial 
magmas. Integration of the settling velocity over tfie time period 
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of the cooling history gives the distance some mineral Is. capable 
of settling. This measure of crystal mobility can be used to 
check whether the Apollo 12 olivine basalts can be genetically 
related by olivine fractionation In a magma body with cooling 
properties commensurate with the observations on 12002. 

Table 2 lists the experimentally determined compositions of 
the residual liquid in 12002 as a function of temperature. Values 
at 50°C intervals were interpolated from graphs of composition 
versus quenching temperature. (Trends In the residual liquid 
compositions were very similar for all experiments cooled at 
<50®C/hour. ) - The densities and viscosities of these residual 
liquids have been calculated according to the models of Bottinga 
and Weill (1970, 1972). Some calculations, however, require 
extrapolation of density and viscosity to lower temperatures 
(<1200°C) than those covered by the models, In these cases the' 
density and viscosity were calculated for the particular compos- 
ition at high temperatures and extrapolated to lower temperatures. 
Density was taken to vary linearly with temperature; the natural 
logarithm of viscosity vms varied linearly with Inverse tempera- 
ture. The density extrapolation Involves relatively small changes 
In value and, therefore, probably small errors. The Ahrrenlus 
extrapolation probably tends to underestimate the viscosity at 
lower temperatures, thus leading to overestimates of settling vel'- 
ooitles. However, as shown below for plcritio liquids, most of 
the settling Is over by the time the body has cooled to 1200°C, 
so these errors are taken to be minimal here . 
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The density variation of the residual liquid Is quite araall, 
decreasing from, 3.29 to 3.26 g/oo hetween 1350°c and 1050"C. 
Evidently the change In chemistry favoring decreased density Is 
nearly compensated by thermal contraction. The variation In. 
viscosity (Figure 3) from ~10 poise near the .llquldus to more 
than 1100 poise near the solidus is quite dramatic by contrast. 
Also shown In Figure 3 as functions of temperature are the degree 
of oystalllnity ($) of 12002 estimated optically and supported by 
crude chemical mass balances-, the maximum size of olivines and th« 
density of the olivines based on average core and rim composltlone 
The settling velocity (V^) .of a sphere In a liquid is given 
by Stokes' Law as 

Y s: (X) 

2 

where g is the acceleration of gravity (1.63 m/seo on the moon), 
V Is the radius of the sphere, Ap Is the density difference 
between the sphere and the liquid, and n is the viscosity of 
the liquid. Using this formula, settling velocities of olivines 
can be calculated as a function of temperature for the data 
shown in Figure 3. The results of these calculations are shown 
In Figure A where the temperatures are put on the time axis In a 
manner appropriate to the cooling history of 12002. The 
settling velocity rapidly increases with decreasing 
temperature from the llquldus to -1250'’C in response to 
the p dependence of velocity on Increasing crystal size. 

At lower temperatures the crystal has nearly attained Its 
maximum size and the effects of logarithmically increasing 
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Viscosity cause a rapid decrease In settling velocity. The 
area under the settling velocity curve corresponds to the 
total distance the grovrlng crystal settles, which In this 
case Is ~16 meters. 

This simple calculation by Stokes' Law is not directly 
applicable to cooling magma bodies which are suspensions of 
non-spherloal crystals. Shaw (1965) and Shaw ^ (1968) 

discuss some of the problems related to settling in silicate 
liquids. Shape factors (K) developed by HcNown and Halalka 
(1950) can be used to correct symmetrical shapes for deviations 
from sphericity. However, for the case of nearly equant 
olivines produced at slower cooling rates no substantial cor- 
rection is needed -(K = 1) . The fact that crystals are settling 
through a suspension of other crystals introduces non-trivlal 
changes which are difficult to treat exactly. The viscosity 
of a pure liquid will be significantly different from that 
of a suspension of particles in the same liquid. Rosooe (1952) 
has extended the Einstein model for the bulk viscosity of a 
auspenslnn to the case of unequal grain- else distribution. 

When this model Is applied to Stokes' Law, equation (i) becomes 

^RE = “ 1.35<t)2-5 (2) 

where Vbj, la the Roscoe-Elnsteln settling velocity, K is the 

AE< 

Shape factor and $ is the degree of crystallinity. When 
equation (2) is Integrated ovei* the cooling time there Is a' 
sharp decrease, in the settling distance (relative to the simple 



Stokes Law Calculation) to -5-7 meters as shown In Figure 
4. The presence of suspended crystals raises the effective 
viscosity and impedes settling. 

An empirical expression for the settling of particles 
in a suspension has been developed by Lewis, Gilliland . and 
Bauer (1949), When this expression is applied to Stokes’ 

Law, equation (1) becomes 

\o, ■ (3) 

where is the Lewls-Gllllland-Bauer settling velocity. 

Equation (3) is probably a more correct model than (2), 
since (2) is based upon the effective bulk viscosity of the 
suspension, rather than the effective viscosity experienced 
by a particle in the suspension. However, when is 
Integrated over time the settling distance Is 4.1 meters, 
which Is not very different from the integration of 

(5,7 meters). Thus we take the distance which olivine 
crystals may settle in a 12002 melt during the deduced 
cooling episode to be approximately 5 meters. 

This estimate of ~5 meters is based upon a cooling 
rate which is Intlally 0.8“C/hour but progressively slows 
down and upon olivines which nucleate at the liquidus and 
reach 2 mm maximum dimensions. If the cooling rate Is assumed 
to be a constant 0,8®C/hour, the settling distance based on 
changes from 5-7 to 5-3 meters. If the initial cooling 
rate is increased to’ 1.5°C/hour, the settling distance 


based on V™ drops to 2.T meters. Obviously the calculation 
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is much more sensitive to changes in the initial cooling 
rate than to changes in the form of the cooling curve 
because most of the settling occurs during the Initial 
cooling of the magma. The calculation Is also' very sensi- 
tive to the size of the olivines. Pew crystals' In the 
natural specimen attain dimensions of 2 mm. Host olivines 
have a maximum dimension less than 1 mm. The smaller 
crystals have a computed settling distance of ~1 meter. 

It Is clear that for the sort of cooling history inferred 
for 12002 the olivine crystals should have sunk out of 
any hand specimen-sised volume of melt in which they 
nucleated, no matter which calculation scheme is used. 

An entirely different conclusion is reached when, the 
calculations, simulate the cooling of the vltrophyre, 12009 
(Donaldson et al., 1975). Here the cooling rate after 
eruption is more uncertain (1,500° C/hour ) . The olivines 
are plaby and skeletal, which requires alteration of the 
shape factor (K' = 2) and density terms (Ap' ” 4p/2) . 

Strong composition gradients In the liquid make it 
difficult to estimate a meaningful value for the 
viscosity. However, even with these uncertainties the 
magnitude of the settling dlstonoes ranges only from 1 to 
3 ram. Obviously these olivines do not settle out from their 
sites of nuclcation and thus aoeoiding to these ealeulatlons 
picritic liquids may certainly be preserved at very fast 
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QooXlng rates. 

Another Indication that these calculations give reasonable 
results In lltaltlng cases can be found by calculating settling 
distances for pyroxenes in 12002. Experimentally grown pyro-’ 
xenes did not attain the dimensions of the natural pyroxenes 
so a representative maximum dimension of 2 nw was taken from 
natural pyroxenes in thin section. An inferred size versus 
temperature plot is shown in Figure 3- Pn** the cooling history 
of 12002, pyroxene settling In 12002 would amount to 3 to 5 ca- 
The results of the settling calculation appear in Figure k. 
Pyroxenes would not be expected to settle far from their sites of 
nucleation. Whatever minor loss did occur from a particular 
volume would be readily compensated for by pyroxenes settling 
into that volume from above.* The great contrast between the 
settling distances predicted for olivines and pyroxenes is a 
result of the very large settling velocities of olivine near 
ISOO^C where viscosity and crystallinity are low. whereas 
pyroxenes begin to grow and sink only after viscosity and 
crystallinity are substantially increased. As a result we 
should expect to find in the 12002 magma body considerable 
variation in the bulk rock chemistry as a result of -olivine 
movements but little with respect to pyroxene. We noted above 


*(Slrailar calculations also indicate that mobility of metal 
spherules, spinel, and plagloolase is very limited as a result of 
the small crystal sizes and large i in the case of plagloolase.) 
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that the major-element chemistry of the olivine basalt suite 
13 described by olivine control without significant pyroxene 
covariation. Apparently, simple olivine settling Is at least 
qualitatively consistent with the differentiation which occurred 
in the natural magma both on the basis of chemistry and the 
experimentally determined physical properties of the magma. 

The previous discussion has served to Illustrate crystal 
settling calculations but must be substantially modified to 
calculate quantitatively differentiation by olivine settling. 

This is because as the crystals move they encounter new 
environments which are cooling more rapidly or more slowly 
than where they started. (The settling distances are 
large relative bo hand specimen size.) We are required to 
consider more than Just the cooling history at one point. To 
do this we must now investigate what can be said of the 
size and cooling history of the magma from which the samples 
came, and the probable sample distribution within that magma 
body. 

MAQMA BODY CHARACTIRXSTICS 

An Important conclusion about the suite of samples can be 
drawn from considering the implications of the strong p.Qsi-tlye 
correlation of normative olivine with plagloolase grain size shown 
in Figure 2. In a magma bods which differentiates by crystal 
settling, olivine accumulates at the bottom and is depleted 
at the top of the cooling unit. Cooling by heat loss to the 
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surroundings is more rapid at the margins and leads to a finer 
grain size at the margins relative to the interior of the body. 
These points are schematically Indicated in Figure 2 for a thin 
magma body, e.g., flow or sill. Notice that only in the bottom 
portion of a coollnp; unit can one find a positive correlation 
of grain size and olivine abundance. The point is well 
Illustrated in terrestrial examples such as the prehistoric 
Makaopuhi lava lake which has been described by Moore and Evans 
(1967, pp. 208-9). A question of some importance which this 
model must then face is: what has happened to the top of the 

magma body in which the Apollo 12 olivine basalts were generated? 
We will show from settling calculations that a substantial por- 
tion of this missing material will be depleted in olivine and 
have the composition of a cotectlo residual liquid. O'Hara 
et al . (1975) have repeatedly pointed out- the quantitative 
importance of cotectic residual liquid composition as a soil 
constituent. Perhaps the soil is the sink for this missing 
material, A possible explanation may be that meteorite impact 
selectively degrades the top portions of thick, differentiated 
plerite flows. Only the basal section (including a lower chill 
margin and some of the olivine cumulates) survives as samples 
large enough to be returned as hand specimens. Whether this 
degradation oooured In the time Interval between flows at the 
Apollo 12 site or In the ~3 b.y. since the last flow there Is 
conjectural. The answer to this question depends to a certain 
extent on what the stratigraphic order of flow types is 
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inferred to be and on the rate of soil production. 

It should be noted that another mechanism is known for 
producing a positive correlation of" olivine abundance with 
grain size. DreVer and Johnston (1967a) have noted this 
relationship in some terrestrial sills where the correlation 
Includes the top as well as the bottom of the cooling unit. 

If this situation applies, rather than the one analogous to 
the prehistoric Makaopuhi lava lake depicted in Figure 2, we 
are spai’ed the embarrassment of finding the missing top of the 
Apollo 12 magma body. However, these cases are clearly related 
to magmas which were emplaced with substantial proportions of 
Intratellurlc phenocrysts since all parts of the body are more 
enriched in olivine than either the top or bottom chill margins, 
unlike the Makaopuhi case where the top of the lake is depleted 
in olivine relative to the chill margins. Bhattacharjl (1967) 
has shown experimentally that Intratellurlc phenocrysts may 
be concentrated In the Interior of flows and sills, relative 
to their chill margins, by flowage differentiation. Furthermore, 
Slmkin (1967) has shown that the flowage differentiation mechan- 
ism is responsible for the distribution of olivine In the 
intrusions such as those Illustrated by Drever and Johnston 
(1967a) to a much greater extent than simple olivine settling. 

An important feature of bodies In which olivine Is distributed 
by flowage differentiation rather than simple settling Is that 
olivine abundance correlate’s with olivine crystal size (in 
addition to the plagioclase grain size which Is not controlled 


EEPRODUCIBILITy OF THE 
OEIG-INAL PAGE IS POOR 
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by flowage differentiation) . 'Shis Is clearly seen in Slrakin's 
{ 1967 ) profiles and Is spectacularly sketched by Orever and 
Johnston U967b)- Olivine grain size variations of this 
magnitude (factor of >2-3) are not seen In the Apollo 12 
piorlte suite and what little variation is present does not 
correlate with olivine abundance (Table 1). This suggests 
that we are not dealing with flowage differentiation in the 
Apollo 12 piorlte suite. Furthermore, it will be shown below 
that settling would substantially modify a homogeneous initial 
distribution of Intratellurlc phenoorysts in an Apollo 12 
piorlte magma so that the final distribution would be quite 
unlike that observed. An Initial distribution of phenoorysts 
concentrated toward the body center such as produced by flowage 
differentiation would make the resemblance even poorer. Sample 
12009 was examined by Drever et al. (1972) for preferred 
orientation In its olivines such as might be expected fro® 
flowage differentiation. This quenched sample' would be the 
most likely place to preserve such effects without later 
reorientation by settling. They did not regard their results 
as indicating any significant effect. We conclude the flowage 
differentiation model probably does not apply to the Apollo 12 ■ 
piorites and will proceed with the simple settling model. 

If this suite is a sampling of the bottom portion of a- 
ooollng unit, two rather Important conditions apply to the 
history of the suite. The heat loss leading to solidification 
was principally through the adjacent magma bottom and, as a 
result of cooling from below, the 3ower portions of the magma 
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body were gravitationally stable with respect to convection 
which would, then, not occur. These two conditions allow one 
to proceed with 8lng>lified thermal calculations~a conductive 
boundary (rather than fixed temperature or radiative surface 
appropriate to a flow top) may be used and convection may be 
neglected to a first approximation. 

We may use the experimentally determined cooling history 
of sample 12002 to calculate how much insulation is necessary 
to hold this cooling element of magma to a certain initial 
cooling rate, i.e., how far away from the basal surface of 
heat loss the magma element must be to have that cooling rate. 
Estimates of the distance from this nearest surface of heat 
loss can then form the basis of an estimate of the dlmanaiona 
of the magma body. Any perturbations Introduced by surfaces 
of heat loss other than the nearest one will increase the 
cooling rate. If there are two surfaces of beat loss com- 
bining to give the observed cooling rate then those surfaces 
must be farther away to be equivalent in cooling power to 
a single surface of heat loss. Therefore, assuming only one 
nearby surface of heat loss (the base) and assigning all the 
conductive heat loss to that surface forms the basis for a 
minimum estimate of magma body size. The single-surface 
approximation will only be valid if the other surfaces of heat 
loss are several times as far away as the nearest one. 
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Therefore an estimate a few times the distance to the nearest 
surface of heat loss can be taken as a minimum estimate of 
the dimensions of the magma body. 

An analytical treatment of this case which considers 
heat of fusion as well as a finite melting Interval la due to 
Jaeger (1957), who Improved upon the treatment of Winkler 
(19^9). The solution for the^ distance from the base (x) 
has the form: 

X = / 5 (H) 

4t'^ 

where £ Is a proportionality constant given by a rather 
complicated expression Involving the temperature range 

of solidification and the ratios of thermal dlffuslvltiea 

« 

of magma and country rock. The dlffusivitles of magma 
and bedrock were taken to be the same (0.01 cm^/seo) except 
that the magma liberates latent heat of fusion (80 cal/gm) 
over a range of solidification from 1325-1000®C which gives 
an effective dlffuslvlty of the magma of only 0.005H cm^/seo 
(vis., Jaeger's equation (1) on p. 307, 1957) for a specific 
heat of 0.29 cal/gm -deg. The bedrock was taken to be 
initially at 0*C. To solve (It) for x, the value of the 


cooling rate (-j^) determined at the llquldus from olivine 

nucleatlon density of ~0.8°C/hour was used (Walker et al. , 

1976a). This value applies to the first 50°C or so of 

cooling (Walker et al., 1976b) and gives a value of £ In 
2 

(It) of -1.7°c«cm /sec when combined with the other constants 
above. The result Is that 12002 Is calculated to have 
crystallized about .90 m from the basal chill margin. 

It is instructive to apply this same approach to the 
vitrophyre 12009 to determine its distance from a cooling 
aurface. Donaldson et al. (1975) noted on the basis of 
crystal morphology and texture that the cooling rate upon 
eruption may have been greater than 100°C/hour. Walker 'et 
al. (1976a) suggested the rate may have been as much as 
500°C/hour. Using these cooling rates and the thermal 
constants used above, adjusted for the lower llquldus 
(1250°C) of 120Q9, one calculates that x is 3-7 cm. Sample 
12009 then Is clearly part of the chill margin of its magma 
body, a conclusion In accord with Its petrography. 

It may be concluded then that the magma body of 12002 must 
have been several times a meter tpick In minimum dimensions 
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according to the considerations discussed above. We point 
out that this conclusion Is not dependent on accepting the 
notion that the suite comes from the bottom portion of a 
cooling unit. It la merely the solution, that comes from a 
conductive boundary condition regardless of configuration. 
Analytical treatments which Involve constant temperature 
heat-loss surfaces are also available (e.g., Carslaw and 
Jaeger, 1959, chapter 11) vfhich would be appropriate to the upper 
surface of a lava flow. We _ find that this treatment places 12002 
about 2.5 meters from such an exposed surface. 12002 was 
calculated to be closer to the conductive surface because 
the insulating country rock in such a configuration produces 
the same cooling rate with less Intervening solidifying magma 
to act as an Insulator. In either situation the 12002 magma 
body must have been at least several meters thick. It could 
be much lai’ger also. ' 

An Important conclusion of our experimental study of 12002 
(Walker et al . , 1976a) was that the cooling rate decreased 
with time during the solidification. The thermal calculations 
Just described also indicate that the cooling rate should 
decrease. Specifically for 12002 the cooling rate Is calculated 
to decrease by a factor of 3 between crystallization of 
llquldus olivine and groundmass and the solidification is 
calculated to take 3-ll weeks for the given initial cooling 
rate. We had previously guessed the rate might have decreased 
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by as much as a factor of 10 during the solidification (which 
might take a month or two) . This was based on extrapolation 
of our experimental observations which were at cooling rates 
too rapid to cover all but the Initial stages of 12002 
crystallization. 

Additional Information on the cooling rates of the other 
samples In the suite would further our understanding of the 
dimensions of the magma body. Unfortunately no direct comparison 
of these samples with controlled cooling experiments Is yet avail- 
able, principally because the coarser samples require longer to 
cool than Is presently feasible experimentally. Extrapolation of 
our experimental olivine nuoleatlon density and pyroxene "swallow- 
tail" curves which Indicate cooling rate is not useful since much 
Of the olivine in the coarser samples Is cumulus and the coarser 
pyroxene morphologies do not have recognizable "swallowtails." 
Fortunately, some progress can be made on the problem by 
considering groundmass plagloclase grain size. 

ORAIN SIZE, COOLING RATE, AND MAGMA BODY SIZE 

Although there are many complications which perturb the relation 
of grain size to cooling history and hence distance from the margin, 
sufficient terrestrial observations have been made to allow some 
cautious statements about the connection of the plagloclase grain 
size variation in our Apollo 12 suite and the size of the magma body. 

Figure 5 shows plots of mineral grain size versus the measured 
distance from the margin of the body for a number of terrestrial 
occurences of different minerals. As expected, grain size 
coarsens away from the body margins. What, then, is the : 
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relationship that we may use to relate grain size to distance 
from the margin? To a first approximation the grain size could 
be considered to be linearly related to distance from the margin; 
however. If we consider the concavity towards the distance axis 
In Figures 5*^, ®, and f to be significant we might wish to 
formulate the relation as follows: 

size a x^ (5) 

where p Is an exponent less than or possibly near'one. It Is' 
possible to estimate p from theoretical considerations as well 
as from Figure 5. From equation (^) we saw that distance (x) 

Is inversely proportional to the square root of the ooollng 
rate. This result Is actually quite general and applies to 
many solutions of the cooling problem other than the specific 
one considered above. The relation of cooling rate and grain 
size is not so straightforward since it Involves a complex 
Interplay of diffusion, nucleatlon, and growth rates. However, 
we found experimentally that the spacing of growth Instabilities 
in pyroxene, for example, displayed a linear relation with 
cooling rate on a log-log plot CWalker et al., 1976a, fig. 10). 
The slope of the line was -i; 

log (size) » -i log (cooling rate) + b (6) 
which when combined with (i|) gives p = | In (5). We do not 
wish to claim that exponent ^ Is general or even correct. It 
Is consistent with the relations In Figure 5 but only marginally. 
If at all, more so than a linear relation. 
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The s’aln size variation of plagioolase In the Apollo 12 
olivine basalt suite may now be used to establish limits on 
the magma body size. Since we know the position of 12002 from 
experimental cooling studies (which give the cooling rate) and 
heat flow calculations (which give the position for that cooling 
rate) we can calibrate equation (5) for the suite; 

2 

plagslze In ram = .004 (distance In cm)'^ (7) 

or 

plagslze In mm = .089 (distance in m) (8) 

depending on choice of exponent. The parallel studies on 12009, 
the vltrophyre, afford a cross check on which exponent Is 
appropriate. Using (7) and data in Table 1, 12009 Is calculated 
to be j cm from the margin; and using (8) the distance Is 2^ cm. 
It would seem that (8) the linear relation, gives slightly 
better Internal consistency. Lane (1903b) suggested theoretical 
reasons why the relation might be linear If equation (4) applies. 
We adopt (8) which gives the most conservative size estimates. 

The coarsest samples represent those furthest from the 

% 

base of the body and they are calculated to be -11 meters from 
the base by (8). [The distance would be -38 meters by (7).] 

The magma body must be at least a few times this distance in 
thickness so that heat loss from the top surface does not 
contribute to a more rapid cooling history than Implied. The 

gafabros come from a magma body that was probably at least 
3_Q .meters In thickness . Although flows of this thickness are 
known (Schaber et al . , 1976) there Is also evidence that lunar 
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lava flows are commonly thinner, Brett (1975) summarized this 
evidence and argued by similar reasoning connecting chemical 
kinetic and heat flow data that most lunar lava flows were 
probably less than 10 meters in thickness. Although no 
kinetic data on these coarse gabbros was available to Brett and 
hence there is no direct conflict, coarse gabbros from Apollo 
15 gave solutions placing them in surprisingly thin flows 
(e.g., 15065). An analysis baaed on plagloolase grain size 
would probably place them In. substantially thicker flows. 

Brett's analysis was based on Zr partitioning between ulvfisplnel 
and llraenlte which Taylor et al. (1973) have shown to be related 
to cooling rate. The temperature range covered by this analysis 
Is large— down to 850®C. Winkler (1949, p. 563) noted that when 
temperature falls substantially below the temperature of intrusion 
(e.g., to 850°C here) cooling rate differences across an Intrusive 
are small. If the Zr partitioning Is responding well down to 
these subsolldus temperatures then the reason coarse and fine- 
grained samples give similar cooling rates by Zr partitioning Is 
understood. In this case, cooling rates determined from features 
developed within the solidification range will be a more reliable 
indicator of position within a magma body than features developed 
In the subsolldus. Therefore v;e tentatively prefer to use the 
analysis based on plagloolase grain size even though It Implies 
substantially thicker magma bodies than previous studies 
indicated. Recognition of uncertainties, especially those 
associated with possible boundary conditions such as layers of 
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soil of low thermal conductivity, may make the results of the 
different calculations more compatible. Table 1 lists 
distances calculated on the basis of equation (8) and these ' 
are used In Figure 8. 

MAGMA DIFFERENTIATION 

We now wish to explain quantitatively the distribution of 
olivine in the Apollo 12 olivine basalt suite on the basis of 
gravltatlve olivine settling Into the bottom portion of the magma 
body which the suite of samples represent. The explanation will 
be satisfactory if it can produce the curve of Figure 8 which 
gives olivine concentration as a function of distance above the 
base. The simplest interpretation of this curve Is that 12009, 
the vitrophyre at the chill margin, represents the magma bulk 
composition In which the accumulation has occured. In this 
interpretation 12002 would be a phenocryst-enrlched liquid 
rather than a liquid composition.* 


*The amount of cumulus olivine suggested from Figure 2 and 
Table 1 is about 9? based on an average normative olivine 
content for 12002 from two different analyses. Our analysis 
of 12002, 57 (the cample crystallized experimentally) suggests 
that it has only more normative olivine than 12OO9. The 
difference In analyses may reflect heterogeneity In olivine 
phenooryst content as noted above for 12004. Furthermore, from 
the data in Table 1 of Walker et al. (1976a) it can be calculated 
that 3 -^% of cumulus olivine In 12002 would explain the difference 
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It will be assumed that 12009 Is a better guide bo initial 
snagma' composition. This assumption and other alternatives Will 
be considered later. 

The calculation of olivine settling profiles Is done by 
computer. A size is postulated for the magma body. At the 
llquldus temperature olivine crystals are allowed to nucleate 
at arbitrary even spaoings. The magma cooling is oonputed 
and the crystals grow and settle. The computation for each 
crystal stops when 'the velocity becomes negligible. The 
olivine conftentration In the final settling profile is computed 
by considering how close the final positions of the crystal 
are relative to their initial spacing. This follows the 
convenient formulation of Pujil (197*)). 

-Settling velocities are computed from equation (3). The 
input parameters as a function of temperature are determined 
from the experimental work on 12002 as shown in Figure 3. 

Olivine rddlus is not taken from Figure 3> which gives the 
maximum experimental olivine sizes. Instead olivine was 
allowed to grow to sizes representative of the olivines 

In composition between the experimental liquidus olivine (FOyg) 
and the olivine actually present in the sample (FOyg). A few 
percent cumulus olivine In 12002 makes little difference in the 
cooling rate cited above since the oumulua olivine can be accounted 
for by a few large crystals. The nucleatlon density, therefore, 
changes little; moreover, the cooling rate estimate la rather in- 
sensitive to nucleatlon density in this range of cooling rate.e. 
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petrographlealiy observed in the suite (Table 1) with a temp- 
erature dependence of size as shown In Figure 3. Various 
sizes were used. Also, for some calculations It was 
assumed that olivine was present In suspenslcn at the time of 
intrusion and settling of these phenoorysts, without further 
growth, was calculated. 
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For calculation purposes the magma body was assumed to be 
a flow with upper surface at 0®C and the heat flow solution of 
Jaeger (I968) was used to determine temperature as a function 
of time and position: 


T = erf 


x+a 

2^/kt 


erf 


3a-x 

2 v1?t 


erf 




2)4ct 


(9) 


where x Is the distance from the midplane of the flow with 
thickness 2a, Tj, is the extrusion temperature ( 1250 “c for a 
magma with I2009 properties), t Is time, erf is the error 
function,. and k is the thermlal dlffusivlty (0.01 omVsec). 
This solution does not consider liberation of latent heat and 
the thermal properties of the magma and bedrock are assumed 
to be the same. We shall see later that certain results of 
the calculations allow use of a solution which does Include 
the effects of liberation of latent heat. 

In order to make computing time manageable, but at the 
same time keep the final computed position of the crystals 
numerically stable, it was necessary to make the time steps 
in the Iterative calculation adjustable. This Is especially 
important when the sinking crystal encounters a regime of 
rapidly changing settling velocity. The procedure adopted 
was to compare the velocity at the start of a time step with 
that at the end of the time step at the new, tentative, po- 
sition of the crystal. If the positions of the crystal com- 
puted with the different velocities disagreed by more than 


500 microns, the crystal was not advanced to the new position 
and the time step was halved until the computed positions at 
the end of the time step differed by less than 500 microns. 

If the difference in final positions was less than 50 microns, 
the time step was doubled for the next step and tested again. 

It was found that decreasing the maximum permissible difference 
to 100 microns did not change the final calculated positions 
of the crystals. 

Figure 6 shows some calculated settling profiles for a 
spectrum of magma sizes and for various different ultimate sizes 
for the growing olivines. Olivine concentration is plotted as 
a function of position in the bodies. One obvious feature of 
these profiles is that larger bodies differentiate more than 
small ones as would be expected since larger bodies cool more 
slowly and allow more settling before solidification. For 
6a and b where crystals grow during solidification, a l-meter 
flow cools too rapidly to experience any differentiation 
whereas a 6 -meter thick flow settles all olivine out of a 
substantial portion of its Interior. For the case of olivine 
crystals growing to a maximum radius of ^ mm (6b) more settling 
occurs than for a maximum radius of j nan (6a). The profile 
is thus very sensitive to crystal size. As we noted before 
In the section on crystal settling, most of the settling 
occurs at temperatures very near the llquldus. If the magma 
Is extruded with a suspension of phenocrysts already developed. 
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the settling velocity Is at a maximum at first extrusion 
since these crystals do not have to grow to sizes large 
enough to settle rapidly. Figure 6c shows profiles com- 
puted for the case of phenocrysts of j mm radius extruded 
with the magma and allowed to settle without further growth. 
There Is obviously a dramatic increase In settling compared to 
6a and b even though the crystals in 6fa grow to larger size. 

The profile Is therefore even more sensitive to crystal size 
near the extrusion temperature than to the final growth size. 

One feature exhibited by 6c and common to all the profiles 
we computed for original suspensions of phenocrysts, is the 
roughly linear slope on the portion of the curve near the 
magma bottom. In some oases the curve Is, actually slightly 
concave towards the concentration axis (e.g., fig. 7a), This 
behavior contrasts with profiles calculated for crystals grow- 
ing during settling which are convex towards the concentration 
axis near the magma bottom although some approach linearity. 

It Is to be noted that the olivine concentration versus distance 
curve for the Apollo 12 suite (figure 8) is distinctly con- 
vex towards the concentration axis. The simplest explana- 
tion Is that this suite represents accumulation of crystals 
which mostly grew after extrusion rather than accumulation of 
a suspension of phenocrysts developed prior to emplacement. 

Another feature of the profiles Is that zones of extreme 
olivine concentration develop complementary to zones completely 
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devoid of olivine, since no provision vfas made In the cal- 
culations to avoid mutual Interference of olivine crystals, 
these ollvlne-enrlched. horizons which develop are unrealisti- 
cally thin and concentrated. An olivine concentration factor 
of 5 Is probably about. the realistic physical limit possible 
and the curves are arbitrarily cut off there. No attempt 
was made to estimate realistic thickness for these ollvlne- 
enrlched zones since none appear to have been sampled (no 
dunltea) and our principal concern here Is the section adja- 
cent to the base and beneath the zone of greatest olivine 
enrichment . 


An Important feature of the profiles In figure 6 is obscured 
by the presentation of many different magma body sizes on a 
scaled distance coordinate. For any particular set of condi- 
tions (e.g., crystal size, extrusion temperature, etc.) the 
olivine concentration curve near the magma bottom Is Independent 
of magma size. The curves for different size bodies overlap 
until the curve for the smaller body either develops a highly 
concentrated zone or falls back to lower concentrations. This 
feature Is observed for every spectrum of magma body sizes 
for every condition we calculated. Two examples are shown in 

L. 

figure 7. This result Is intuitively plausible when one con- 
siders that the settling which occurs near the margin will be 
controlled by the heat loss at that margin, For regimes near 
the margin this Initial heat loss Is Independent of the size of 
the magma body. We have argued that the Apollo 12 olivine basalt 
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suite does sample such a regime near the magma base. If this is 
true then we can calculate the settling profile for the lower 
portion of a magma body large enough to Include the Apollo 
12 olivine basalt suite as this lower zone with an analytical 
heat flow solution which includes the effect of latent heat . 

We used such a. solution (Jaeger^ 1957) in the discussion 
above relating cooling rate to distance from the base, and 
it can be applied here also: 

T =■ erf x no) 

2/kF 

*^L ^^^50 C) and Tg (lOOO^C) are the llquidus and solidus 
temperatures; x is the distance from the base; t is the time; 

It is the effective diffusivlty of the magJiia when a latent heat 
of 80 cal/gm is liberated uniformly over the solidification 
range (here taken as 0.0048 om^/sec); and y is a complicated 
expression Involving the solidification temperatures and ratio 
of effective diffusivlties of magma and country rook (with a 
numerical value here of 0.546), 

Settling profiles for the bottom ’of a magma were then 
calculated with this formula fbr a large magma thickness 
(arbitrarily taken as 64 meters)-. These are shown in figure 0(a and b) 
The consideration of latent heat makes a noticeable difference 
to the settling profile. Figures 6a and b show dashed durves 
calculated this way for comparison to calculations without 
latent heat. Settling and concentration at the base is enhanced 
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in these models as compared to those which do not take latent 
heat into account. The inclusion of latent heat undoubtedly 
gives a more realistic result. We repeated the set of 
calculations for a magma with 12002 physical properties 
(iji^a!]^325°C; k=.0054 om^/seo; Y“-456) to check the relative 
merits' of the two magma compositions. The higher llquldus 
temperature of 12002 with slightly more normative olivine 
results in increased settling for any particular crystal 
size. (See curve 8c.) 

Also plotted in Figure 8 is the Apollo 12 olivine basalt 
suite. It la clear that the settling profile for olivines 
groMlng to a radius of | mm In a magma with 12009 extrusion 
temperature (1250“C) and composition gives a very good 
description' of the Apollo 12 olivine basalt suite foi> the upper 
part of the curve. However, the suite seems to have concen- 
trated slightly more olivine In its basal members than would 
be expected for a 12009 magma composition. On the other 
hand, a 12002 magma composition will concentrate too much 
olivine in its basal section to be consistent with the Apollo 
12 olivine basalt suite . This may Indicate either that the 
actual magma was intermediate between 12002 and 12009, that a 
small fraction of the olivine In a 12009 magma was present at 
the time of emplacement as phenoorysts, or that some combination 
of these effects may have operated. In either case, the 
deviations of the suite from the calculated profiles are 


rather small. 
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We conclude that simple olivine settling can give a 
quantitative. Internally consistent explanation for the 
chemical and petrographic properties or the Apollo 12 
olivine basalt suite. 
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TABLE 1 

APOLLO 12 OLIVINE BASALT DATA 



12009 

1^075 

12004 

12002 

12Q2Q 

12018 

12040 

■13035 

s modal olivine 

13.6 

18.2 

11.9 

19.6 

14.3 

14.1 

24.6 

28.7 

t range of Fo 
In olivine 

76-49 

77-51 

76-54 

77-50 

77-32 

76-49 

67-44 

63-50 

t normative 
olivine 
content 

12.3 

17.5 

17.9 

21.6 

20.5 

25.4 

27.4 

30 .,9 

* mineral size 
in mm: 









plagloclase 

,002 

.05 

.06 

,08 

.11 

,16 

.46 

.98 

olivine 

.56 

-77 

.53 

.72 

.60 

.63 

.72 

.51 

4= distance In 
meters above 
base 

.02 

.56 

.67 

.90 

1.2 

1.8 

5.2 

11.0 


* Average values from 1 to 3 sections of each rock. Samples used were; 
12009 , l4j 12075,24; 12004,51; 12002,7,9, and 160; 12020,8 and 15; 
12018,6.77, and 795 1204o,44 and 45; 12035,21 and 22. 

t Average of values from analyses reported by: Compston et al. (1971), 

Cuttltta ^ 8^. (1971), Engel et al. (1971), Grove et al. TT973) , 
Kushlro and Haramura (1971), Maxwell and Wllk (19717T Schoon D971), 
and Willis ^ a]^. (197I) . Table. 10 of James and Wright (1972) gives 
a convenient compilation. 


4 Calculated by equation (8) in the text from plaglooiase size. 



TABLE 2 


'IS 

FRACTIONATED RESIDUAL LIQUIDS 


T°C 

1350 

1300 

1250 

1200 

1150 

1100 

10 SO 

Si02 

43.7 

43.6 

44.3 

45.0 

46.0 

46.7 

46.5 

TiOg 

2.5 

3,0 

3.3 

3.5 , 

3.8 

4.2 

4, 5 

AI2O3 

8.0 

8.5 

9.5 

10-0 

11. S 

12.5 

14.0 

Cr2C>3 

0.70 

0, 85 

l.O 

0.6 

0,3 

0.15 

0.07 

FeO , 

21.9 

21-3 

- 21.1 

20.4 

19,3 

18.9 

19.5 

Mgo 

13.5 

11.4 

9i8 

8.2 

6.1 

4.1 

2.7 

CaO 

8.0 

8,7 

9.5 

10-1 

10. 8 

11,1 

10.8 

KjO 

0.05 

0.06 

0.07 

0.07 

0.08 

0.09 

0.10 

NajO 

0.23 

0,27 

.0.30 

0.33 

0. 36 

0.40 

0.45 


98.58 

97.88 

9B.87 

98.20 

98.24 

98,14 

98-72 


Residual liquid compositions interpolated from experimental 
comparisons of residual liquids produced at cooling rates 
< 50“c/hour. K 2 O and Na20 values necessary for the viscosity 
and density calculations are estimated from initial bulk con- 
tent and degree of crystallinity. This disequilibrium frac- 
tionation trend may be compared to the equilibrium sequence 
In table 1 of Walker et al. (1976a)- 






velocities shown in Figure 
FlKure 4 

Calculated settling velocity for olivine and pyroxene ■ 
is shown throughout the 12002 crystallization history as 
calculated by various methods discussed in the text . Note 
that the velocity scale for pyroxene is much expanded relative 
to that for ollvlnev Areas under the curves are a measure 
of how far olivine or 'Pyroxene can settle. These distances 
are given in parentheses after the symbol. Indicating which 
formula was used to calculate the curve. 

Figure 5 

Plots of linear dimensions of minerals versus their 
distance from cooling surfaces. Literature sources are: 

Lane (1903a) j Moore and Evans (1967.) > Queneau (1902), and 
F. Walker (19*10). Queneau 's data is plotted as the square 
root of area reported. F. Walker's data for plagloclase 
In the olivine diabase horizon of the Palisades sill is 
indicated in separate symbols since K. R. Walker (1969) 
has shown that this material (and what is above) is involved 
In a second magma injection which perturbs a simple cooling 
history. Queneau evidently avoided this material in his 
analysis. 

Figure parts a, b, and c Include data only from the 
marginal facies of the solidified bodies. The plots are 
more nearly linear than parts d, e, and f where the observations 
include material at the center of the body. The marginal 


portions of d, e, and f can be considered approximately linear. 
The curvature of plots d, e, and f reflects the additive 
cooling effects of both margins on the grain size In the body 
center. 

Parts a, o, and d include data from more than one section 
across the bodies and indicate the sort of lateral variations 
thay may be found in a body. All the minerals compiled here 
are thought to have crystallized after magma emplacement (with 
the possible exception of those in part f) and the patterns 
here should nob be Influenced by flowage differentiation 
during magma emplacement . 


Figure 6 

Calculated profiles of olivine settling for flows of 
various thickness. Vertical position coordinate is sealed 
for total thickness which Is indicated in meters by the 
circled numbers. The olivine concentration scale on the 
abscissa is relative; unity is the Initial value for each 
of the magmas which are assumed to be homogeneous liquids 
or homogeneous crystal-liquid mixtures at time of emplacement. 
Equation (9) gives temperature history. Dashed curves in 
lower portions of a and b give results from equation (10) 
for comparison of effect of latent heat on a 6*1 meter thick 
flow. Results in a and b are for crystals which grow during 
settling as contrasted to o where phenocrysts of ^ mm radius 


ire present at time of emplacement ^and do not grovf during 
settling. 
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SOLID SOLUTION AND PHASE EQUILIBRIA ON 
THE JOIN ANORTHITE-SILICA 

J . Longhi 

J. P. Hays (both at: Dept, of Geological- 

Sciences, Harvard Univ. , Cambridge, 

Mass. 02138) 

Melting experiments have been conducted 
x^ith divitrified glass compositions along 
the Join anorthit e-silica. Microprobe 
analysis of the run products shows good 
agreement of liquid compositions with 
those determined by Schairer and Bowen 
{Bull. Soc. Geol. de Finlande, 19^7). 
However, microprobe analysis of the crys- 
talline phases shows significant solid 
solution in both anorthite and silica 
phases. At l485'’C there is approximately 
2 wt ? excess silica in anorthite coex- 
isting with liquid; at 1422° C the solid 
solution increases to 5 wt ^ and reaches 
a maximum of approximately 8 wt % at the 
eutectic, 1368° C. There is also approx- 
imately 5 wt ? solid solution of calcium 
and aluminum calculated as anorthite 
present in eristobalite at the eutectic. 
Silica-rich mixtures of anorthite glass 
and crystalline tridymlte were equili- 
brated to determine the equilibrium 
stability relations of .the silica poly- 
morphs. Inversions of tridymibe to 
eristobalite v/ere obtained in less than 
a week at 1374° and 1333°C. Thus tridy- 
mite has no liquidus stability field 
along the Join and very limited solid 
solution (<1 wt ^). We conclude that 
preferred substitution of impurities and 
not temperature is the dominant factor 
controlling the presence of tridymite 
versus eristobalite. On the basis of, 
derivative structures (Beurger, Am. 
Miner., 1954) we’-can predict- that ’ cris- * 
toballte will be the stable phase in Ca 
and Na-rlch igneous environments, whereas 
tridymite will be stable in K— rich 
igneous environments. 
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FELDSPAR FLOTATION AND LUNAR CRUST FORMATION 

WALKER, D. and HAYS, J. F., Hoffman Laboratory, 

Harvard University, Cambridge, Mass. 02138 

The lunar crust is thought to have formed as the result 
of feldspar flotation during an early melting event 
involving a substantial volume of the moon. An attempt 
was made to show that feldspar would indeed float during 
such an event. Large anorthite crystals ( 'v.S mm) were 
placed beneath a silicate glass of Mg/Mg+Fe =0.83 {*) 
representative of liquid in which plagioclase accum— 
ulation^is thought to have occurred. In less than 3 hburs 
at 1300 G, the crystals rose to the top in a Pt cruci- 
ble 3 cm deep equilibrated in air and in a Mo crucible 
1.5 cm deep equilibrated in an H^/CO^ gas stream of - 
log PO 2 = -10,9 (below Fe/FeO) . '^Previous attempts 

to observe feldspar flotation were probably hampered 
by use of the smaller crystals characteristically grown 
from experimental charges. ' These results suggest that 
lunar crustal formation by feldspar flotation is possible 
without special recourse to differential sinking of 
plagioclase versus mafic minerals or selective elutriation 
of plagioclase. 

(*) SiO^ 47.5; TiO^ 0.9; Cr202 0.2; AI 2 O 2 20.7; 

"FeO" 4.5; MgO 12.0; CaO 12.3; K 2 O 0.4; Na20 0.7 
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